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FOREWORD 


This Program User's Guide discusses the Nuclear Thermal Propulsion (NTP) engine 
system design features and capabilities modeled in the Nuclear Engine System Simulation (NESS): 
Version 2.0 program (referred to as NESS throughout the remainder of this document), as well as 
its operation. NESS has been upgraded to include many new modeling capabilities not available in 
the original version delivered to the NASA Lewis Research Center in December 1991, see 
Ref. 1-0. NESS's new features include: 

• An improved input format 

• An advanced solid-core NERVA-type reactor system model (ENABLER II) 

• A bleed-cycle engine system option 

• An axial-turbopump design option 

• An automated pump-out turbopump assembly sizing option 

• An off-design gas generator engine cycle design option 

• Updated hydrogen properties 

• An improved output format 

• Personal computer operation capability 

Sample design cases are presented in this user's guide that demonstrate many of the new 
features associated with this upgraded version of NESS, as well as design modeling features 
associated with the original version of NESS, discussed in Ref. 1-0. 
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CONSULTATION 


Please direct all questions or comments associated with the Nuclear Engine System 
Simulation (NESS) program pertaining to the overall model to: 

Mr. Dennis G. Pelaccio 

Science Applications International Corporation 

21151 Western Avenue 
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Indianapolis, IN 46250 

(317) 576-5697 

FAX: (317) 841-7599 


Issues related to the engine system reactor model should be addressed to: 

Dr. Lyman J. Petrosky 
Westinghouse Electric Corporation 
Advanced Energy Systems 
Waltz Mill Site 
Madison, PA 15663 
(412) 722-5110 
FAX: (412)722-5234 


General model operation questions and program distribution issues should be addressed to: 

Mr. James Walton 
NASA Lewis Research Center 
Nuclear Project Office 
21000 Brookpark Road 
Cleveland, OH 44135 
(216) 977-7047 
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1.0 INTRODUCTION 


An accurate, standalone, preliminary Nuclear Thermal Propulsion (NTP) engine system 
destgn analysts tool is required to support cutrent and future Space Exploration Initiative (SEI) 
propuls, on and vehtcle design studies. Currently available NTP engine design models are those 
eveloped dunng the NERVA program in the 1960s and early 1970s and are highly unique to that 
destgn (see Ref. 1 - 1 ) or are modif, cations of current liquid propulsion system design models. To 
ate, TP engine-based liquid design models lack integrated design of key NTP engine design 
eatures, such as tn the areas of reactor, shielding, multipropellanl capability, and multiredundam 
pump feed fuel systems. Add.tionally, since the SEI effort is in the initial development stage a 
robust, verified NTP analysis design tool could be of great use to the community. 

™ s effon de ''= l °!* d ™ accurate, versatile NTP engine system design analysis program 
( ), known as the Nuclear Engtne System Simulation (NESS) program, to support ongoing and 

future engtne system and stage design study efforts. In this effort, Science Applications 
In emanonal Corpora, ion s (SAIC) NTP version of the Expanded Liquid Engine Simulation 
( S) program was modified extensively to include Wesdnghouse Electric Corporation's near- 
term and next generation solid-core reactor design models, ENABLER I and ENABLER II 
reactor designs, respectively. The ELES program has extensive capability to conduct prelimimuy 
system design analysis of liquid rocket systems and vehicles. The ptogram is modular in nature 
and is versatile ,n terms of modeling state-of-the-an component and system opdons as discussed in 
Rets. 1-2 and 1-3. The Wesdnghouse reaclor design model, which were integrated in the NESS 

program, are based on the near-term and upgraded version of the solid-core ENABLER NTP 
reactor design concept, see Ref. 1-4. 

Hus program is now capable of accurately modeling (characterizing) a complete near-term 
or next generation solid-core NTP engine system in great derail, for a number of design opdons, in 
an efficient manner. The following discussion summarizes die overall analysis methodology key 
assumptions, and capabilities associated with die NESS, presents example problems and 

compares the results to a related NTP engine system design. 


1-1 



2.0 ENGINE SYSTEM MODEL 


Th 1S section discusses the overall NTP engine system design and performance prediction 
methodology and the unique model input options associated with NESS. To better understand the 
operation of NESS, it is important that the operator be familiar with the ELES program which is 
discussed in detail in Refs. 1 -2 and 2- 1 . 

2.1 Overall Analysis Method 

The NESS flow logic is essentially the same as the ELES logic detailed in the ELES 
Programmer's Manual, Ref. 1-3. A simple summary of the analysis procedure is shown in Figure 
2-1, and a detailed flow chart is g.ven in Figure 2-2. Many portions of the code are iterated two or 
more t.mes to improve accuracy. The key inputs include the thrust level, FVAC, reactor type, 
IREACTR, and engine cycle type, KCYCLE=1 for gas generator (GG), =3 for expander, or =7 
for bleed cycle. Also important are the chamber pressure and temperature, PC and TCHAMBER, 
respectively, flow paths (bypass fractions NFF and BYPTUR), nozzle configuration, NOZTYP 

and KOOLNZ, turbopump type, 1PTYPE, reactor scaling factor, FALPHA, and the number of 
propellant feed legs, NTP A. 

Once an input file has been formulated and read in by NESS, the first step is to initialize 
propellant properties from the libraries of propellant data stored in the code. These properties will 
be recalculated at many different code locations and for many different conditions throughout code 
execution. The ideal performance is initially estimated based on known chamber pressure and 
temperature, and nozzle area ratio; the boundary layer and divergence efficiencies are calculated at 
this time and an estimated delivered specific impulse (Isp) is found. This estimate is used to 
calculate a reactor flowrate. The nozzle heat load is estimated as 1% of total reactor power, and 
this heat load, Isp, and flowrate are passed to the reactor design portion of the code, ENABLER, 
for calculation of reactor fuel and overall operating characteristics. The generic NESS ENABLER 

reactor design module can be configured to represent either an ENABLER I or ENABLER II 
NTP reactor design, see Section 3.0. 

The reactor inlet pressure and temperature are now used to calculate the cycle pressure 
schedule. During the pressure calculations, the nozzle barrier cooling requirement is also 
calculated along with the regen cooling requirements. Now that all engine efficiencies are known, 
the actual dehvered Isp and flowrate are calculated. The actual nozzle heat load is compared with 
the original estimate and if they are not within 10%. the code loops back to the reactor design 
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portion of the code and repeats all steps up to the point this comparison is made. If the nozzle heat 
loads are reasonably matched but the reactor design has only been performed once, the code loops 
back to the reactor design with the newly calculated Isp and flowrate to improve accuracy. 



Figure 2-1. NESS Program Overview 
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Figure 2-2. NESS Program Flow Logic 
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After the reactor design, performance, and pressure schedule have been completed 
satisfactorily, the code now calculates all cycle flowrates. Tankage volumes, pressures, 
temperatures, and pressurization requirements are calculated next. The temperature schedule is 
determined, and the turbomachinery can now be analyzed. The turbopump assembly (TPA) 
portion of the code calculates the size and performance of the pumps and turbines, and checks for 
cycle balance by comparing pump required horsepower to the turbine delivered horsepower, if not 
balanced, a new turbine pressure ratio is calculated and the TPA design process is repeated. 


Once the TPA design has been completed, the flowrate and temperature schedules are 
recalculated to improve accuracy. Next, component weight calculations for the reactor, 
turbomachinery, nozzle, and all miscellaneous parts (lines, valves, etc) are performed. Mass 
flowrates are calculated one more time, overall engine dimensions are found, and finally, output 
summaries are printed out. When the pump-out (double run) option is selected (see Section 
2.3.1), the entire design process is completed for an engine at reduced thrust level and then a 
second iteration of the entire design at full thrust level is performed beginning with the reactor 
module using some of the values calculated in the first pass (TPA parameters and some weights). 

Flow path schematics of the representative NTP expander, gas generator, and bleed engine 
cycle systems are shown in Figure 2-3. The representative NTP engine systems shown in Figure 
2-3, incorporate dual propellant feed systems in all cases, and boost pumps for the expander and 
bleed cycles only. The representative GG system does not include boost pumps. 


2.2 Major Code Components 

Table 2-1 lists the major code modules along with key flags and input variables. Each of 
these modules is discussed in further detail in the sections following, including both overall 
discussion of the module and how to determine the inputs required. 

2.2.1 Engine Performance 

Engine performance calculations begin with an ideal one-dimensional equilibrium (ODE) 
performance value that is later degraded with loss multipliers. The ideal values for specific impulse 
(lsp) and C star (C*) are calculated by the ODE module of the Two-Dimensional Kinetic Reference 
Program (TDK), Ref. 2-1 , as a function of chamber pressure, temperature, and nozzle area ratio. 
Tables of hydrogen performance data are stored in the subroutine HYDROGEN along with the 
curve-fit equations used to calculate ideal C*. which is a function of temperature and pressure 
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Figure 2-3. Representative NTP Expander, Gas Generator, and Bleed Engine System Cycle Row Paths 






Table 2-1. Key NESS Input Flags and Variables 


Module 

Variable 

Value 

Result 

General Input: 




Cycle Type 

KCYCLE 

= 1 

Gas Generator Cycle 


= 3 

Expander Cycle 



= 7 

Bleed Cycle 

Thrust Level 

FVAC 

- 

Set thrust level 

Chamber Pressure 

PC 

- 

Set pressure 

Chamber Temperature 

TCHAMBER 

- 

Set temperature 

Reactor Type 

IREACTR 

= l 

Use ENABLER 1 


= 2 

Use ENABLER D 

Choose double-run? 

IDBLRUN 

0,1 

If =1, double -run used 

Iterate pump design? 

ITRATE 

0,1 

If=l, iterate pump design 

Bleed cycle solver 

ISOLVE 

= 0 

Input bleed flow fractions 


= 1 

Input turbine inlet temp. 

Input bum time? 

IUSRBRN 

0,1 

If =1, input bum time 

User-defined TPA? 

ISTSET 

0,1 

If =1, input TPA values 

Nozzle: 




Exit area ratio 

EPS 

- 

Set exit area ratio 

Use extension? 

KEXNOZ 

0.1 

If =1, use extension 

Use 3 -portion nozzle? 

NOCTYP 

0.1 

If =1 , use 3-portion nozzle 

Attach area ratio 1 

EPSATT 

- 

Set ext. attach area ratio 

Attach area ratio 2 

EPSAT2 

- 

Set 2nd ext attach area ratio 

Nozzle cooling 

KOOLNZ 

= 2 

Regen cooling of nozzle ext. 


= 3 

Trans-regen cooling 



= 4 

Radiation cooling 



= 5 

Film cooling (GG cycle only) 

Regen Cooling: 



Set turbine bypass flow 

Turbine Bypass Ratio 

BYPTUR 

- 


= 1.0 

Tie tube flow drives turbine 

Barrier Temp. Fraction 

DIFTBF 

- 

Set barrier temperature 

Reactor: 



Original flow paths 

Reactor flow paths 

CONFIG 

= 1 


= 2 

Tie tube flow drives turbine 

Fuel Type 

FTYPE 

= 1 

Graphite fuel 


= 2 

Composite fuel 



= 3 

Carbide fuel 

Support Pattern 

SPAT 

= 2:1 
= 3:1 

Set support pattern 

Nozzle Flow Percent 

NFF 

= 6:1 

Set nozzle/tie tube flows 

Fuel Scaling Factor 

FALPHA 

- 

Set fuel scaling factor 

Tankage: 



Tandem tankage 

Tank Type 

NCTNK 

= 0 


= 1 

Non-conventional tankage 

Pressurization Method 

KGASFL 

= 0 

Cold helium or solid GG 



= 1 

Autogenous 

Turbopump Assembly: 

JCNFIG 

= 1 
= 2 

Gearbox 

Single shaft TPA 

Pump Configuration 



= 3 

Twin TPA in series 



= 4 

Twin TPA in parallel 



= 5 

Multiple feed leg TPA 

Turbopump Type 

IPTYPE 

= 0 
= 1 

Centrifugal Pumps 
Axial Pumps 

Use Boost Pumps? 

JBPFL 

0,1 

If =1, use fuel boost pump 

Number of Feed Leas 

NTPA 

- 

Set number of TPA feed legs 
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only. An ideal Isp at desired conditions is interpolated from these tables. To run the code with a 
propellant other than hydrogen, ODE (or a similar code) must be run to generate the tables of Isp 
data and the C* equations. This data would be put into a new subroutine that is called by the rest 
of the code when appropriate. 

The loss multipliers used to degrade the ideal performance are calculated using standard 
JANNAF procedures. Ref. 2-2, or Aerojet-derived methods. Ref. 1-2. It is assumed that the 
reactor itself has no losses, and therefore engine efficiency is determined by nozzle-related factors. 
The efficiencies ( or losses) calculated by NESS are the nozzle boundary layer efficiency, 
divergence efficiency, and nozzle barrier cooling efficiency. The gas generator bleed efficiency is 
calculated when applicable. A thorough explanation of these efficiencies is given in the ELES 
Technical Information Manual, Ref. 1-2, and the key equations are summarized below. 

The boundary layer loss equation was developed by Aerojet as a result of their experience 
in defining this loss. The equation is as follows: 

ETABL = 0.997 - (In(EPS)/100)*[l.-0.065*ln(0.01*Pc*Fvac)+ 

0.00 1 *(ln(0.0 1 *Pc*Fvac))2] 

where 

EPS = Nozzle Exit Area Ratio 
Pc = Chamber Pressure (psia) 

Fvac = Vacuum Delivered Thrust (Ibf) 

This equation is accurate for engines with a radiation or film cooled nozzle, but does not 
take into account the energy returned to the core flow by a regen-cooled nozzle. In this case, the 
energy lost by the nozzle is retained by the regen coolant flow and fed back into the engine, and 
therefore should not be considered a true loss. A nozzle that is completely regen cooled should 
have a boundary layer efficiency of 1.0, while a partially regen-cooled nozzle, as is typically used, 
should have an ETABL less than 1.0, but higher than that predicted by the above equation. To 
provide accurate modeling of the regen-cooled nozzle option, an input adjustment factor, ADJBL, 
is applied to the efficiency calculated by the above equation. The adjustment factor is applied as: 

ETABL = 1.0 - (1.0 - ETABL)*ADJBL 
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The current value used for ADJBL of 0.2 (code default = 1.0) was determined by 
comparison with Rocketdyne performance values, see Ref. 2-3, which were calculated in much 
greater detail than is possible with NESS. 

The divergence loss is a function of nozzle shape and was derived as curve-fits of the 
information presented in Appendix A of the CPIA document No. 178, see Ref. 2-4. The equations 
are as follows: 

For conical nozzles: 

ETADIV = 0.5 + cos(alpha)/2. alpha = half angle in deg. 

For RAO nozzles: 

ETADIV = 1.0-(l.-C)*[(1.75-RATMLR)/0.75)] 1 7 for RATMLR <= 1.75 
or 

ETADIV = 1.0 for RATMLR > 1.75 


where 

C = constant = 0.945 + 0.01*ln(EPS) for EPS <= 20 

= 0.958 + 0.00566*ln(EPS) for EPS > 20 

EPS = Nozzle Area Ratio 

RATMLR = ratio of nozzle length to the length of a minimum 
length RAO nozzle; an input 

The divergence efficiency can also be adjusted, if desired, with the input factor ADJDIV used as: 
ETADIV = 1.0 - (1.0 - ETADIV)* ADJDIV 

The barrier cooling loss is a function of the amount of coolant fluid needed to maintain the 
nozzle wall temperature below the maximum allowable for the material used. Aerojet chose a 
simplified barrier cooling loss routine consisting of a stream tube analysis which flow-averages the 
performance of the core stream tube with that of the barrier stream tube. The procedure for 
calculating stream tube flow areas and flow rates is detailed in the ELES Technical Manual, Ref. 
1-2. The maximum barrier temperature is input as described in section 2.2.2, and is used to 
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calculate barrier Isp and C*, and ultimately barrier mass flux. The fraction of fuel used for barrier 
film cooling (FFFC) is calculated as: 

FFFC = barrier flowrate/(barrier flowrate + core flowrate) 

The barrier loss (ETABAR) is set at 0.95 and is put into the comprehensive barrier cooling 
loss equation: 

ETAMRD = [(Isp*mdot) core + (Isp*mdot*ETABAR) baiT j er ]/(Isp*mdot) toto i 
where all Isps are ideal. 

This efficiency can be adjusted by the input ADJMRD in the same form as that used for the 
boundary layer and divergence losses. Note that the "barrier cooling loss" is referred to as the 
"mixture ratio maldistribution loss" in the ELES manuals. 

For gas generator cycles, the gas generator bleed efficiency is calculated as a function of the 
bleed nozzle flowrate, pressure, and area ratio. It can be adjusted with ADJGGB in the form: 

ETAGGB = ETAGGB * ADJGGB 

All other efficiencies described in the ELES Technical Manual, Ref. 1-2, were set equal to 
1.0 because of their inapplicability to the nuclear engine; for example, injector or fuel and ox 
mixing efficiencies. 

2.2.2 Nozzle Cooling 

The nozzle can be cooled by a number of methods. The converging portion of the nozzle, 
including the throat, is automatically regen cooled. It is of milled slot construction to upstream area 
ratio of 4 with an adapter of regen tubes connecting the nozzle to the reactor. The remainder of the 
nozzle is cooled by regen tubes, radiation, a cold film of turbine exhaust (GG cycles only), or by a 
combination of these. A detailed explanation of regen cooling calculations is given in the ELES 
Technical Information Manual, Ref. 1-2, and Section 2.2.3 of this report gives nozzle modeling 
options. 

The nozzle regen cooling requirements are based on the nozzle wall material properties, 
chamber temperature, regen coolant flowrate, regen inlet temperature and pressure, and regen 
channel size. The maximum wall material temperature in input as TGWNOM and is the 
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temperature above which the material will begin to degrade. For copper, a common converging 
nozzle material, this max temperature is 1460°R. The 1460°R temperature limit is typical of that 
used for the maximum design nozzle wall temperature for the Space Shuttle Main Engine (SSME) 
which is made of NARLOY-Z, a copper alloy, Ref. 2-5. For the high chamber temperatures 
typical of nuclear reactors, the regen coolant is unable to maintain this maximum wall temperature 
if the fluid on the other side of the wall is at chamber temperature. Therefore, a small amount of 
cool fluid from the regen outlet is dumped into the chamber at the top of the converging nozzle and 
is used to form a cool barrier between the wall and the hot core fluid. The loss in efficiency due to 
this barrier cooling is detailed in Section 2.2. 1 of this report and in the ELES Technical Manual, 
Ref. 1-2. The greater the temperature mismatch between the barrier fluid and the core fluid, the 
larger the cooling loss, and therefore the highest possible barrier fluid temperature should be 
chosen that can still maintain the required material wall temperature. The barrier temperature is 
input as a relation between the core temperature and max wall temperature, TGWNOM. The input 
variable DIFTBF is used as follows: 

T bamer = TGWNOM + DIFTBF*(T core - TGWNOM) 

Ideally, DIFTBF = 1 .0 and the barrier temperature equals the core temperature to minimize 
flow losses. If DIFTBF = 0.0, the barrier temperature is set equal to the max wall temperature. 
For a copper wall with max temperature 1460°R and a core temperature of 4860°R (2700°K), the 
maximum barrier temperature that could still maintain the required wall temperature is 1630_R, 
which means the input DIFTBF = 0.05. A good value for DIFTBF can really only be determined 
by past experience and trial and error; the larger the difference between the maximum wall 
temperature and the core temperature, the lower the value for DIFTBF will have to be. 

Other key regen cooling inputs include the gas wall material thermal conductivity and 
minimum gauge. The land width (WLTHR) and channel width (WTHR) of the regen cooling 
channels at the throat are also important inputs because they will strongly affect the regen pressure 
drop, i.e., small channels => high velocity => large delta P. There is also an option for user-input 
regen pressure and temperature drops, initiated with the flag INDPDT set equal to 1 and DELTAT 
and DELTAP input 

2.2.3 Nozzle Modeling Options 

The user has a number of different nozzle modeling options. The most basic option is to 
set the nozzle extension flag KEXNOZ to zero and have regen slots all the way out to the exit area 
ratio EPS. This type of nozzle is almost never used in practice because of excess weight, and 
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therefore a nozzle extension option is allowed. If the nozzle type flag NOZTYP is set to zero and 
KEXNOZ = 1, an extension will be added to the regen slots. This section extends from area ratio 
EPSATT to EPS, and can be regen, radiation, or film cooled (GG cycles only), with cooling 
option selected with the variable KOOLNZ. The new and final option is for NOZTYP=l, which 
models a three-section nozzle made up of regen slots, regen tubes, and a radiation cooled 
extension. The user must set KEXNOZ = 1, KOOLNZ = 2 (regen tubes in portion 2), and area 
ratios EPS, EPSATT (attach point of second section) and EPSAT2 (attach point of third section). 
Figure 2-4 shows the three nozzle modeling options and key input variables. 

The regen slot portion of the nozzle extends out to an upstream area ratio of 4 where it 
attaches to a nozzle/reactor adapter that is made of aluminum regen tubes covered by a l oad -bearing 
casing of aluminum. The weight of this assembly is calculated in the reactor weight subroutine, 
and is included in the reactor pressure vessel weight. 

Material density and strength are input for the converging nozzle, first nozzle extension, 
and second nozzle extension with RHCSTR, RHONZE, RHONZ2 and SIGCHM, SIGNZE, 
SIGNZ2, respectively. The minimum thicknesses of the two possible extensions are input as 
TNZMIN and TNZMN2. The volume of material used for the regen slots is calculated and the 
total converging nozzle weight is a function of this volume, the density of the material used for 
each region of the slots, and total surface area. The weight of the regen tubes is a function of the 
maximum pressure in the tubes, surface area, and material density, strength, and minimum gauge. 
The radiation-cooled extension weight is simply a function of surface area and material density and 
thickness. 

2.2.4 Reactor 

A solid core, ENABLER-type reactor design module was developed by Westinghouse 
Electric Corporation and integrated with ELES to form NESS. The reactor design is made up of 
two segments: the first calculates fuel requirements and reactor operating conditions, the second 
calculates approximately 30 reactor component weights along with key reactor dimensions. NESS 
provides hydrogen data, Isp, core flowrate, and nozzle heat load to the reactor module 
(ENABLER) for its calculations. In return, ENABLER provides the reactor inlet and tie tube oudet 
conditions needed for pressure and temperature schedule analysis. A detailed discussion of the 
reactor model can be found in Section 3.0. 


2-11 



NOZTYP - 0: Singi«-S«ct>on Nozzfa 


No CxtMion - 
ftagon Sloes oat 



KOOLNZ • 2 


0tmntm»4 

Note: Thu option is Mttom uMd boesuoo at ** one mm * watghc. 


NOZTYP - a Standard Noxzla 


Raqmi, NmHhIov Trmo^opon 



AmRMI -4 


“I ‘ 

IPSATT 



Figure 2-4. Nozzle Modeling Options and Key Input Variables 
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_ . Th ' E 1 NABLER reactor design module consisls of *wo disimct reactor modeling options. 

' sl lsENABL£R '■ a near '«™ solid core reactor design based on the NERVA reactor. The 
second option, ENABLER II, provides a more advanced reactor design, with different flow paths 
and scaled fuel, reflect, ng state-of-the-art technology. This option yields reactor designs with 
tgher power dens.ttes and lower weights. By properly utilizing the inputs to the NESS reactor 
destgn module (ENABLER), these NIP reactor options can be represented. 

One key reactor input is the nozzle flow fraction, NFF, which determines the percentage of 
ow going to the tie tubes and to regen cooling. The user also selects the fuel type as either 

TTlr 0 ? ° r Carb ' de USi "® ' he VaiiM ' FTYPE T* r ««or temperature is input as 
CHAMBER, and ts used extensively in the reactor design process, along with determining the 
overall engine performance. 

The input variable SPAT is used to select the ratio of fuel to support elements in the reactor 

and may be tnput as a support pattern of 2:1, 3:1. or 6:1. The support pattern wiU affect the reactor 

we, ght, overall engtne cycle performance, and reactor criticality. The support pattern should 

always be se, as htgh as possible (6:1) to reduce weight. However, for small engines and scaled 

uel engines, this ratio wtll often need to be reduced in onler to achieve reactor criticality Also 

for some htgh pressure expander cycles, the tie tube flow will no, contain enough energy to drive 

the turbtnes and the support pattern can be reduced to increase the tie tube oude, temperature and 
thus add more energy to the flow. 

Two reactor flow path options are currently available for the expander cycle. Standard 
operation utilizes tie tube flow plus some input percentage of reflector outflow to drive the 

“ renector n ° w percentage is input as turbine bypass ratio, BYPTTJR, (see Section 
2.2. 10 for a detailed explanation), and is used to detetmine wha, fraction of the reflector outflow 
ypasses e turbine. The second flow path option is always used with the ENABLER II reactor 
esiSh option and results in dm turbines being driven by the tie tube flow only, with all reflector 

xr° ,he corc -. msnow ^ cm ** ^ * ° ne * *** 

BYPTIIR tn ? rSCt h™ 10 ei, ber choice automatically sets the other variable, i.e., if 
= !' “ ' ,h ' inpu, file ' C0NFIG is automatically se, B 2 by the code and vice versa. 

c^le cTJtfrr '1 ' “ d refleCt0r PrcSSUre " “ d ,hOTfore a lighter-weight pressure vessel. This 
yc e can be difficult to power balance, especially for high chamber pressures, and die typical 

or support pattern of 6:1 must occasionally be reduced to 3:1 in order to increase the tie tube 

emperature and add enough energy to the turbine drive fluid for successful cycle balance. 
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This support ratio reduction yields a much higher support mass that often outweighs the reduction 
in pressure vessel mass. 

An important reactor input is the fuel scaling factor, FALPHA. It provides a means of 
including state-of-the-art technology in the reactor design. If FALPHA is equal to 1.0, the 
resulting reactor design represents an ENABLER I NTP reactor configuration which is typical of 
NERVA technology. A value for FALPHA less than 1.0 simulates the advances made in fuel 
element design and corresponds to an ENABLER II NTP reactor system design. This multiplier 
is applied to all fuel element cross-sectional dimensions. To scale the fuel element length as well, 
the variable I Ft - must be set to a value other than zero so that it overrides the input core length LC. 
This length scaling is used in the form: 

if LEL >0: LC = FALPHA * LEL 

The power per fuel element is input with PMW, and is actually the power per a 52 inch length that 
is reduced based on the scaled element length. The minimum value physically possible for 
FALPHA is 0.5; in practice, however, the recommended minimum value is 0.67. 

The flag IREACTR is used to determine the reactor option to be used for the run, i.e., 
ENABLER I or ENABLER II. If IREACTR = 1, the ENABLER I method is used, and the fuel 
scaling factor (FALPHA) is automatically set to 1.0. The variables CONFIG and BYPTUR are 
input by the user to set the flow paths. For IREACTR = 2, the user inputs FALPHA (default = 
0.67) and the flow paths are automatically set to the second option, CONFIG = 2 and BYPTUR = 
1.0, so that only the tie tube flow drives the turbines. 

During reactor design, a criticality check involving the number of ZrH loaded tie tubes is 
made, and a warning printed out if criticality cannot be achieved. For small engines 
(< approximately 35,000 lbf) this warning will typically be issued, but the user should further 
evaluate the design. The criticality determination is based on a comparison of the minimum fuel 
volume required for each support pattern with the actual calculated volumes. Some small engines 
may not pass this fuel volume check, but could actually achieve criticality if other parts of the 
reactor, such as the reflector, were designed differently. The detailed design analysis that would 
be required for an alternate reflector design, for example, is not available in ENABLER. 


An estimate of fuel life is calculated based on the allowable mass loss rate. For scaled 
fuel, the support to fuel element ratio (SPAT) will often need to be reduced in order to achieve 
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criticality. This ratio reduction will yield a smaller minimum fuel volume requirement that can be 
now be met with the scaled fuel volume. 

As can be seen in the worksheet, the user can input a number of variables related to heat 
pickup in various sections of the reactor, as well as several fuel element characteristics. The 
channel coating thickness at various locations may be input with the variables ZRCI, ZRCO, and 

ZRCH. The pressure vessel material properties such as specific gravity and allowable stress may 
also be input. 


The fraction of possible ZrH loading in the tie tubes is input with the variable FZRH and 
varies with engine size. This fraction has not been studied extensively, and therefore the values 
presented in this section are merely suggested values. For small engines with a diameter less than 
approximately 25 inches, FZRH will typically be equal to 1.0. For large engines with unsealed 
fuel (diameter > approximately 34 in.), FZRH should be 0.0. For scaled fuel engines, compare 
the core volume with the standard NERVA 35 in. x 52 in. length reactor core, for which 

FZRH=0.0, and scale FZRH accordingly. A core diameter of 31 in. uses an FZRH of 
approximately 0.4. 

The weight of the control safety rods required by the reactor during launch is calculated by 
ENABLER, but is not added into the nominal reactor weight. During standard operation, these 
safety rods are placed in the reactor during launch and are then discarded upon achieving a safe 
orbit. In both the reactor summary and final engine summary, the safety rod weight is listed after 
the reactor or engine weights, and a total launch weight is then calculated. 

2.2.5 Auxiliary Components 


The category 'auxiliary components" consists of instrumentation, a pneumatic supply 
system, thrust structure, gimbal, and reactor cooldown assembly. Previously in ELES, some of 
these component weights were calculated as a percentage of the total engine weight, some were a 
function of thrust only, and some were not calculated at all. Also, these weights were originally 
calculated assuming a standard liquid rocket engine rather than a nuclear rocket engine. 

A report issued by TRW, Ref. 1-1, includes equations for various nuclear rocket engine 
auxiliary component weights. These correlations relating component weight to reactor power were 
developed as curve fits of NERVA-type reactor data. The TRW equations applicable to the 
ENABLER-typc rocket engine design have been programmed into NESS and which arc: 
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Instrumentation: weight = 166.9 + 0.00743*P - 1.64E-7*P^ 

Pneumatic Supply System: weight = 751.6 - 0.00208*P + 2.35E-6*P^ 

Reactor Cooldown Assembly: weight = 238.1 + 0.0254*P - 8.04E-7*P^ 

Upper Thrust Structure: weight = 786.25 - 0. 1868*P + 5.2E-5*P^ 

Lower Thrust Structure: weight = 492.9 + 0.091 1*P + 1.463E-6*p2 

where P = power in MW 

The upper and lower thrust structures are combined into the "thrust mount" weight. The other 
three weights make up the "support hardware weight". 

Although these equations provide a useful starting point for auxiliary component weight 
calculations, they represent NERVA-era technology rather than state-of-the-art designs. To 
account for advances in technology, weight multipliers arc input that decrease these weights to 
values more in line with current engine designs. The instrumentation multiplier, CXWINST, is 
left at 1.0. The pneumatic supply system weight was compared with similar system weights on 
current engines, such as the SSME, and was found to be extremely high, see Refs. 2-5, and 2-6. 
It should be noted that the TRW pneumatic supply system weight correlations assume that the 
complete pneumatic supply is part of the NTP engine system, while for the SSME the main supply 
is located in the Space Shuttle. This is one major contributor to the weight difference as well as the 
higher pressure and lighter weight components associated with today's systems. Therefore, the 
pneumatic system weight multiplier, CXWPNEU, is input as 0.25. The reactor cooldown 
assembly multiplier, CXWTNKAS, is input as 0.9 to account for technology advances. The thrust 
structure multiplier, CXWTHM, is set to 0.9 to allow for lighter weight materials and improved 
technology. If NERVA-era technology is desired, all above multipliers should be input as 1.0. 


2.2.6 Materials of Construction 

The NESS user is allowed to select the material of construction of all major subsystem 
components. Standard library tank materials include 6061-T6 aluminum and 6A1-4V titanium, or 
the user may input density, strength, and minimum gauge for a previously undefined material. A 
discussion/comparison of candidate cryogenic tank materials is given in the ELES Technical 
Information Manual, Ref. 1-2. The input worksheet includes a table of the most common engine 
materials along with their densities and strengths. This data is typically used for valves, nozzles, 
lines, and regen channels, and the user may input data for any unlisted material desired. The 
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nozzle designs also require input of minimum material thicknesses. The turbine blade strength and 

density, as well as an overall TPA density that is used in pump and turbine weight calculations, can 
also be input. 

2.2.7 Tankage 

The mam tankage options in NESS are either tandem tankage, in which fuel and oxidizer 
are stacked on top of each other to fit within a common shroud, or non-conventional tankage, 
where the user selects the number of tanks as well as their shapes and placement on the stage. The 
tandem tanks option should probably not be used for nuclear thermal rockets because they use only 
hydrogen as propellant, and may carry only a very small amount of oxidizer for use with a gas 
generator. The tandem tank model automatically calculates an oxidizer tank weight even if the 
amount of oxidizer carried is very small or zero, and this tank is sized to fit in the tank shroud with 
a diameter based on the size of the large fuel tank. The non-conventional tankage design option 
should give a better estimate of actual tank sizes. 

The tank sizes for both tank geometries are dependent on amount of burned propellant, 
ullage fractions, acquisition system design, residual propellant, propellant boiloff, and autogenous 
pressurization. The approach taken in sizing the propellant tanks is as follows: 

a. 

b. 

c. 

d. 

e. 

f. 


These tank volumes are now used to determine pressurization requirements and update initial 
estimates. 

The large variety of possible tandem tank configurations is shown in the ELES Technical 
Information Manual, Ref. 1-2, along with the equations used to calculate many of the tank 


Amount of fuel burned is input; calculate amount of oxidizer burned in GG if 
necessary. 

Add weight of autogenous pressurization requirements to each prope llan t 

Calculate the tank free volumes using the propellant densities and ullage fractions 

Calculate propellant residuals and acquisition device volumetric displacement based on 
tank free volume estimate 

Calculate tank surface area as needed for heat transfer calculations to determine 
propellant boiloff 

Total tank volume is now calculated as the sum of the above volumes: burnt propellant 
ullage, residuals, boiloff, autogenous pressurant, and acquisition devices 
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dimensions and volumes. All tanks can be cylindrical, spherical, or elliptical (CSE tanks), and the 
non-conventional tankage option allows toroidal tanks as well. Non-conventional tank weights are 
calculated from an ideal tank weight through the use of a tank non-optimum factor, which is 
defined as the ratio of actual tank mass to ideal tank mass. The ideal tank mass is based on tank 
wall thickness and size. The actual mass includes any additional material required for weld lands 
and fittings. For conventional tanks that require feedlines, supports, pressurization, and a 
propellant management device, a tank non-optimum factor of 1.7 is suggested. Different factors 
are recommended for different tank types, and these factors arc listed in Table 7.3. 1. 1 in the ELES 
Technical Manual, Ref. 1-2. The tank nonoptimum factor is input as the variable CXWTNK. 

When preparing inputs for tankage design, the user must first set the variable NCTNK 
equal to either 0 for tandem tanks or 1 for non-conventional tanks. If tandem tanks are chosen, the 
user now determines such factors as arrangement of propellant (fuel forward or aft, etc), common 
or separate dome tanks, monocoque or suspended arrangement, tank head ellipse ratio, tank dome 
orientation, safety factor (SFFLTK, SFOXTK, SFPRTK), and tank material (MTNKFL, 
MTNKOX, MATPT). 

To use the non-conventional tank option, the user should first sketch the arrangement of 
tanks and engines on the stage. The total number of non-conventional tanks is input with 
NTANKS (includes ox, fuel, and pressurant), up to the maximum of 15 tanks. The type of fluid 
contained within each tank is input with the variable INTNK1, where an input of 1 is for ox tanks, 
2 is fuel, 3 is pressurant. For example, if two ox and two fuel tanks are desired, input INTNK1 = 
1,2,1, 2. This indicates that tanks 1 and 3 are ox tanks, and tanks 2 and 4 are fuel tanks; retain this 
same numbering scheme when defining the remaining tank parameters. Input the tank ellipse ratio 
for each tank with ELTNK1. The tank type is selected as either CSE or torus with the variable 
KTANK1. The angular location of each tank gives its relative position on the stage and is input as 
TANGLl. Tank radial location indicates the tank distance from the center of the stage, RADLOl = 
4*1.0 places all four tanks at the edge of the stage and RADLOl = 0 places a tank at the center of 
the stage. Engine angular and radial locations are input similarly with the variables ENGAN 1 and 
ENGRD1, with a maximum of five engines allowed. The material for each tank is selected with 
the variable MATNK1. Tank safety factors are input with SFTNK1, and tank weight multipliers 
are input with CXNCT1. More input variables for each tank geometry are contained in the 
worksheet, see Appendix A. 

The forward and aft skirt length inputs are actually input as fractions of tank lengths. For 
tandem tankage, both aft and forward skirt lengths should be input as 1.0 to form a skirt fully 
covering both tanks. To shroud non-conventional tankage, the forward skirt should be set to 0.0 
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and the aft skirt length should be 1.0. This will yield a skin that covers all tankage and is as long 
as the tallest non-conventional tank. DMOTOR is used to input the stage diameter. 

2.2.7. 1 Tank Heat Transfer. For the long duration missions proposed for nuclear rockets, 
tank heat transfer and insulation arc important aspects of vehicle design. A detailed discussion of 
this area is provided in the ELES Technical Information Manual, Ref. 1-2, and includes 
information on optimizing insulation thicknesses. 


NESS offers four possible tank heat transfer scenarios: ignore tank heat transfer, external 
boundary exposed to conductive source, worst case solar radiation, and ground hold ice formation. 
The desired option is selected with the variable KHXOPT. The most common options are either to 
ignore heat transfer (when tank design is not important) or worst case solar radiation. The solar 
radiation option requires input of insulation characteristics, space hold time, flight time, average 
orbital distance from earth, and earth and solar heat flux parameters. The insulation is typically 
composed of a layer of spray-on foam insulation (SOFI) plus a multi-layer insulation (MLI) 
blanket. The density, thermal conductivity, and thickness of each type can be input. Table 2-2 
lists these values for a variety of types of MLI. 


Table 2-2. Multi-Layer Insulation Data Comparison 


MLI 

Configuration 



K * 

(m 3f ) 


Watts 

(m-K) 

BTU 

(hrft°R) 

DAM/DBL silk net 

19.7 


45.2 

2.82 

4.5xl0' 5 

2.5x10-5 

DAM/Tissue glass 

39.4 


51.9 

3.24 

2.5xl0 5 

1.4x10-5 

SAM Crinkled 

15.7 

40.0 

14.6 

0.91 

4.7xl0- 5 

2.6x10"5 

DAM/SGL Nylon Net 

31.5 

80.0 

53.8 

3.36 

3.0x10-5 

1.7x10-5 

DAM/Dexi glass 

23.6 

60.0 

58.8 

3.67 

5.0x10-5 

2.8x10-5 

DAM Crinkled/Tissue glass 

23.6 

60.0 

31.1 

1.94 

7.0x10-5 

3.9x10-5 

Superfloc 

11.8 

30.0 

13.8 

0.86 

4.5x10 s 

2.5x10 5 


2.2.7.2 Propellant Tank Pressurization. Propellant tanks can be pressurized by cold 
helium gas, a solid gas generator, or autogenously. The method of pressurization is selected with 
the variables KGAS, KGASFL, and KGASOX as shown in the worksheet. The selection of a 
propellant acquisition device, either some sort of bladder or surface tension device, has a strong 
effect on the pressurization calculations. An extremely detailed discussion of tank pressurization is 
presented in the ELES Technical Information Manual, Ref. 1-2. 
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When cold gas pressurization is selected, KGASFL.KGASOX = 0 and KGAS = 2, the 
user also inputs the cold helium storage pressure as PICG and the helium tank final pressure 
fraction, FPULCG, where a value less than 1.0 indicates a blowdown tank. If KGAS is set equal 
to 1 instead of 2, a solid gas generator will be used which requires fairly extensive user inputs 
regarding solid fuel characteristics and bum rates (see worksheet). If KGASFL, KGASOX are set 
to 1, the tanks will be pressurized autogenously. This option has an advantage over helium 
pressurization when the additional weight of the evaporated propellants is less than that of the 
helium storage vessel, as occurs in pump fed stages with low NPSH requirements. The propellant 
used in autogenous pressurization of the hydrogen tank will be bled off from the turbine exhaust 
for all engine cycles. Because only a small amount of oxidizer is used in the GG cycle, the 
oxidizer tank is assumed to be pressurized with cold gas to reduce cycle complexity. If autogenous 
pressurization is selected, the pressurizing oxygen flow will be bled off from the oxidizer pump 
outlet flow. 

2.2.8 Propellant Pressure/Temperature/Flowrate Schedules 

The propellant pressure, temperature, and flowrate are calculated at key points within each 
engine cycle. The pressure schedule is calculated "backwards", beginning with the chamber 
pressure and working back up through the cycle using input and calculated pressure changes. The 
temperature and flowrate schedules begin at the tank outlet and flow down through the cycle to the 
reactor inlet conditions. NESS can handle expander, gas generator and bleed cycles. 

For all engine cycles, the tank outflow is divided into tie tube and regen/reflector flow 
based on the input flow fraction, NFF. The regen flow is used to cool both the nozzle and 
reflector, with a small amount bled off to form a cool barrier inside the nozzle. 

As can be seen in the expander cycle flow paths shown in the schematic in Figure 2-3, the 
reflector outflow can be either dumped directly into the core or used to run the turbine. Reflector 
outflow going to the turbine is mixed with the tie tube flow, and turbine inlet temperature is 
calcu lated by an energy balance of tie tube and reflector flows, i.e. 

^turbine inlet = ((T*mdot*Cp) re fl ector + (T*mdot*Cp)Q C tubel/ 

(Cp*mdot) tur bine inlet 
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where 


T = temperature 
mdot = mass flow rate 

Cp — specific heat coefficient for constant pressure 

Turbine outflow is dumped into the reactor core, with a small amount bled off for autogenous 
pressurization if needed. 

The key pressure assignments for the expander cycle are the turbine and reflector outlet 
pressures. The reactor inlet pressure and temperature are calculated by the reactor model, and are 
therefore known. The tie tube pressure drop is fixed at 250 psid, and the reflector pressure drop is 
25 psid. These pressure drops are typical of solid-core reactor systems, based on past 
Westinghouse NTP reactor design experience. Pressure drops could be higher for large, high- 
heat-load NTP reactor designs. The reflector, turbine, and reactor pressures are related by the 
following list which includes the key pressure variable names and descriptions, along with some 
key pressure cycle assumptions: 

PTURBI, PTURBO = turbine inlet and outlet pressure, respectively 
PREFI, PREFO = reflector inlet and outlet pressure, respectively 
PTTI, PTTO = tie tube inlet and outlet pressures, respectively 
PREGI, PREGO = regen inlet and outlet pressures, respectively 
PCI = core inlet pressure 

PVLVFO = main valve outlet pressure 
TURBPR = turbine pressure ratio 
DELTAP = regen pressure drop 

PTURBO = PTURBI/TURBPR 
PREFO = PREFI - 25 
PTTO = PTTI-250 


For the expander cycle with partial or no turbine bypass (some or all reflector flow goes to 

the turbine), the reflector (PREFO) and tie tube outlet (PTTO) pressures are set equal to the turbine 

inlet pressure, PTURBI. The turbine outlet pressure is set equal to the reactor inlet pressure 
PTURBO=PCI. 
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Once the reflector outlet pressure is known, the reflector inlet pressure, which equals the 
regen outlet pressure, can be calculated so that the regen cooling analysis can be performed and all 
other pressures in the cycle can be calculated. For multiple feed leg TPA designs, the ind.v.dual 
turbine flow rates are multiplied by the number of legs to accurately calculate the pressures. 

For all cycles, the main valve outlet pressure is normally calculated as the reflector outlet 
pressure plus the pressure drop across the regen and reflector, but the valve pressure must be high 
enough to allow for all pressure drops across the tie tubes and turbine. Therefore, the valve outlet 
pressure is set equal to the maximum of the required tie tube inlet pressure and the reflector outlet 
pressure plus regen and reflector pressure drops, i.e., 

PVL VFO=MAX IMUM((PTTOt-250) ,(PREGO+DELTAP)) 

Another option for the expander cycle is to set the input variable BYPTUR equal to 1. 
which sends all reflector flow direcdy into the reactor so that the turbines are driven by the he tube 
flow only. The user must exercise caution when choosing this flow option as the be tube flow will 
occasionally not have enough energy to power the turbines, especially at high chamber pressures 
and when the reactor support pattern is set at 6:1. A support ratio of 3:1 or 2:1 will yield htgher he 
tube temperatures, and therefore more energy to drive the turbine, but the turbine tnlet 
will also be increased, and may exceed the accepted temperature limits of approximately 1400 R. 
A lower support rado will also substantially increase the weight of the reactor. 

When this oprion is selected, the reflector (PREFO) and turbine outlet (PTURBO) 
pressures are set equal to the previously determined reactor inlet pressure (PCI). Tbe tie tube oudet 
pressure is set equal to the required turbine inlet pressure, and the de tube inlet pressure allows for 
the fixed 250 psi pressure drop across the tie tubes. The valve outlet pressure is calculated as 
shown above, and once the valve outlet pressure is set, the pump discharge pressure can be 

determined. 

The gas generator bleed cycle flow schemadc shown in Figure 2-3 uses small amounts of 
oxidizer and fuel to feed the gas generator that drives the tutbine. The turbine exhaust is ether 
dumped overboard through a small bleed nozzle or is dumped into the main nozzle for film 
cooling. Although this exhaust dump results in a performance loss, the GG cycle has the 
advantages of reladvely simple cycle design (TPA and regen design are not coupled) and lower 
pump discharge pressures. Since the turbine is powered by the GG, the reflector and tie tu 
flows are dumped directly into the reactor core. PREFO and PTTO are set equal to PCI, and the 
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remaining calculations proceed as usual. The tie tube inlet pressure, PTTI, is now calculated and 
compared with the valve outlet pressure, PVLVFO, and if PTTI is greater than the valve outlet, 
PVLVFO is set equal to PTTI. This adjustment will occur whenever the regen pressure drop (delta 
P) is less than the fixed tie tube pressure drop of 250 psid. As in the expander case, once the 

reflector oudet pressure is known, the regen cooling analysis can then be performed and all other 
pressures calculated. 

The bleed cycle is analyzed using the same pressure assignments as those used for the GG 
cycle. For the bleed cycle, a small amount of flow from the cold, high pressure propellant pump 
outlet flow is tapped off to combine with the hot, lower pressure flow bled off from the reactor 
chamber exit region to drive the turbine. This cycle is analyzed using the same pressure 
assignments as those used for the GG cycle. The chamber bleed flow undergoes a pressure drop 
as it travels through the line to the mixer. The cold bleed line pressure at the mixer inlet is set equal 
to the chamber bleed line pressure at the mixer inlet to prevent flow backup. The cold bleed flow 
also undergoes a pressure drop as it travels through the lines. These line pressure drops are 
determined by the inputs CPLINH, CPLINC. The remaining pressure drop in the cold bleed flow 
required to match the hot bleed pressure occurs in the cold bleed valve and is calculated 
automatically. After the two flows are mixed, further pressure drops occur across the turbine inlet 
line and turbine throttling valve; these drops are determined by the fractional inputs CPLINT, 
CPVLVT. The remaining pressure schedule for the bleed cycle is calculated by the same methods 
used for the GG cycle. The temperature at the mixer outlet is calculated using an energy balance: 

(Cp*T*mdot) mixer outlet = ( c P*T*rndc>t) hot biced + ( c P*T' : " , «dot) co | d Weed 


To evaluate the bleed cycle, the user must select from the two solver options using the input 
variable ISOLVE. If ISOLVE equals 1, the user inputs the turbine inlet temperature with 
TURBTIN and the code determines the mass flow fractions of hot and cold bleed flow required to 
provide that temperature. The other option is for the user to set ISOLVE equal to zero, input the 
hot and cold mass flow fractions FRACHB and FRACCB, and have the code determine the turbine 
inlet temperature. In practice, the first method will be selected most often because it eliminates the 
extra step required by the second method, namely evaluation of the output to determine whether the 
calculated turbine inlet temperature falls within acceptable limits, and if not, another run must be 
made with mass flow fractions adjusted appropriately. 

For all engine cycles, tank outflow is equal to the core flowrate plus the nozzle barrier 
flowrate, autogenous pressurant flowrate, and gas generator or bleed flow. 
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2.2.9 Propellant Properties 


Propellant properties are required over a very wide range for the variety of models used in 
NESS, including both gas and liquid phases. The approach used to obtain these values is to begin 
with a known value of the propellant property at some reference point, and then scale that value to 
some other condition based on empirical or theoretical correlations. The exceptions to this method 
include hydrogen and helium, which require separate, extensive data bases from which desired 
values are interpolated. A detailed discussion of the methods used to determine property data can 
be found in the ELES Technical Information Manual, Ref. 1-2. Hydrogen data is stored in the 
routine H2DATA. 

A computer program was recently developed at NASA Lewis Research Center to provide 
parahydrogen thermal and transport properties that match the National Bureau of Standards (NBS) 
parahydrogen data, see Ref. 2-7. The NBS data represents the most recent compilation of 
hydrogen properties available in the nation. The program NBSPH2 has been incorporated in the 
NESS program. The routine PH2 was developed to match NBS data exactly across the pressure 
range 29 to 2320 psia (0.2 to 16 MPa) and temperature range 24.8 to 54,000°R (13.8 to 3,000°K). 
The routine includes data tables for density, thermal conductivity, viscosity, Prandtl number, speed 
of sound, enthalpy, and specific heat. All data is stored in SI units and a routine was Written to 
convert the data into the English engineering units required by NESS. 

NESS will occasionally require hydrogen data outside of the pressure or temperature range 
available in NBSPH2. In this case, the original hydrogen data routine, H2DATA, will be called 
instead. For pressures in the range from 2320 psia to 2600 psia, the old routine H2DATA will be 
called with a pressure of 2600 psia and the new routine, NBSPH2, is called with a pressure of 
2320 psia. A linear interpolation using the actual pressure is then performed to find average 
property values. This interpolation was added to prevent fluid property discontinuities in the above 
pressure range. The original (non-NBS) hydrogen data is also used for fluid at the high 
temperatures associated with carbide fuel reactors. If more extensive NBS hydrogen data becomes 
available, that data could be incorporated into the NESS properties at a later date. 

An option exists in ELES that allows for user-defined propellants, which requires that the 
user input certain propellant properties and then select a propellant from the existing ELES library 
that the new propellant is most similar to. The code next evaluates this new propellant performance 
based on comparison with the chosen similar propellant. This option is set up for use by non- 
nuclear, chemical bipropellant propulsion systems, and therefore cannot be used for reactor 
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designs without major code modification. Hydrogen is currently the only propellant with full 
performance data tables programmed into the code, and the current method of determining Isp is 

different than that used for bipropellants and may not be compatible with the ELES user-defined 
propellant evaluation method. 


2.2.10 Turbopump Assembly 


The purpose of ,he "rtopump assembly (TPA) model is lo determine the size, weight, and 
performance of all pumps and turbines for expander, gas generator, and bleed cycles. The code 

can evaluate both centrifugal and axial lurbopuntps. NESS offers the following turbomachinety 
configurations: 3 


L comlSfr dr ' Ving “ 8earb ° X WhiCh P ° WerS an oxidizer “<* Pump on a 

2. Single turbine driving ox and fuel pumps on a common shaft. 

3. Twin TPAs, series drive fluid flow. 

4. Twin TPAs, parallel drive fluid flow. 

5. Multiple propellant feed leg TPA - each leg is identical and sees 1/NTPA of the flow 


The desired option is indicated with the input variable JCNF1G. If the multiple feed leg option is 
chosen (JCNFIG=5), the number of feed legs is input as NTPA. Boost pumps may be included in 

the propellant circuit by setting JBPFL, JBPOX= I . with the boos, pump fraction of total propellant 
head rise input as BPFRFL, BPFROX. 


NESS checks the necessity for pump or turbine staging, allowing up lo four stages for 
cenmfugal pumps, twenty stages for axial pumps, and two stage turbines. To avoid unrealistic 
designs, the code checks the maximum allowable lip speeds and the turbine blade root stresses 
Pump head coefficients and pump and turbine efficiencies are calculaied ftom tables included in the 
program. A parital admission turbine is designed if blade heigh, falls below 0.3 in. The equations 

used to design the centrifugal pumps and turbines are given in the ELES Technical Information 
Manual, Ref. 1-2. 


For htgh flow rates, the low fluid density of hydrogen leads to a high volumetric flow rate, 
a regime for which the multi-stage axial pump is well suited. The axial pump is attractive for such 
applications compared to the centrifugal pump in ierms of wetght, construction, and performance. 
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For this reason, an axial pump option has been added to NESS. A typical axial turbopump 
schematic is shown in Figure 2-5. Axial pump design is selected with the input variable 
IPTYPE=1. Code modifications assumed that an axial pump will not be used for oxygen flow 
(very poor design selection due to the high density of oxygen) and will therefore not be used for 
the gas generator cycle. The logic embedded in NESS to design an axial pump is displayed in 

Figure 2-6. 

The performance calculation methods for the axial pumps are essentially the same as those 
for the centrifugal pumps. Key axial pump design modeling considerations are that the maximum 
number of stages allowed is twenty, and the specific speed (SS) at which the pump will stage is 
3200 (vs. 800 for centrifugal pumps). The pump head coefficient is interpolated from data tables 
containing values based on existing axial pump designs (see Figure 2-7 and Ref. 2-8). The best-fit 
equation used to calculate the head coefficient as a function of main pump specific speed is: 


Head coef. = 0.88237 - 2.3145E-4*SS + 2.3161E-8*SS 2 - 7.7028E-13*SS 3 



Figure 2-5. Typical Axial Turbopump Design 
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Calculate main pump RPM based on NPSP, 
suction specific speed, volumetric flowrate. 

>1 

J"" Calculate total pump head rise. | 

i 

Calculate inducer head rise based on input 
inducer/boost pump fraction 

>1 

Recalculate pump RPM using the new, high inlet 
pressure. 

I 

Design main pump => number of stages, 
diameter, efficiency, horsepower. 


Set inducer diameter and RPM to match that of the 
main pump. 

i 

Calculate inducer performance => specific speed, 
efficiency, horsepower. 

i 

j Total horsepower = main pump + inducer. 

i 

Overall efficiency = weighted average of pump + 
inducer efficiencies, using horsepower as the 
weighting factor. 

i 

Calculate weights and lengths; total pump weight 
I — main pump + inducer. Weight calculations use 
a length to diameter ratio based on existing axial 
pump designs. 


Figure 2-6. Axial Pump Design Logic 
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Specific Speed 

Figure 2-7. Axial Pump Head Coefficient as a Function of Specific Speed 

The maximum allowable tip speed is 1500 ft/sec for hydrogen, above which the pump will stage. 
The axial pump inducer is modeled as a single stage boost pump, with the boost pump flag, 
JBPFL, initialized automatically within the code when the axial pump option is selected. The 
inducer is forced to operate at the same speed (RPMs) and have the same diameter as the main 
pump. Its pressure-head is determined by the input fraction BPFRFL. 

The inducer efficiency is calculated from the existing boost pump efficiency curves in 
ELES, and the main pump efficiency is interpolated from the data shown in Figure 2-8, which is 
also based on existing axial pump design data, see Ref. 2-8. The best-fit equation used to calculate 
main pump efficiency as a function of main pump specific speed is: 



Specific Speed 

Figure 2-8. Axial Main Pump Efficiency as a Function of Specific Speed 
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Main Pump Efficiency = 0.54854 + 1.3501E-4*SS - 1.7544E-8*SS 2 + 5.901E-13*SS 3 
Tlie overall axial pump efficiency is calculated as a weighted average of the inducer and main pump 
efficiencies, using pump horsepower as the weighting quantity, shown as: 

Overall Efficiency = l(HP*efO induccr +(HP*efO pump ]/total HP 

Inducer weight is calculated using the standard boost pump method. The axial main pump 
differs only slightly from the centrifugal main pump weight and is as follows: 

Main Pump Weight = rho*(pi/4)*D 3 *(lVD)*N*f m 


where: 

rho = pump material density (lb/in 3 ) 

D = pump tip diameter (in.) 

(L/D) = pump length to diameter ratio per stage 

N = number of pump stages 

f m = pump material fraction = (0. 12*D + 0.9)/D 

The length to diameter ratios (IVD) per stage for both the main pump and inducer are calculated by 
correlations of the data on length to diameter ratios of existing axial pump designs, see Ref. 2-8, is 
shown in Figures 2-9 and 2-10, respectively. The points on the graph indicate existing design 
values, while the curve defines the correlation used in the L/D calculation as given by: 

OVDlinduccr = 1.992 - 0.23348*D + 0.0106*D 2 

(L/D) main pump = (0.52415 - 0.027 14*D + 0.001 1387*D 2 )*N 

Total axial pump weight combines the main pump and inducer weights. Comparison with existing 
axial pump weights, as shown in Section 2.2.1 1, indicates that a multiplying factor may be 
necessary to bring the pump weight to within the accepted range. 
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Figure 2-9. Axial Main Pump Length to Diameter Ratio (Per Stage) 
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Figure 2-10. Inducer Length to Diameter Ratio (Per Stage) 
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Turbopump lengths are determined so that a feed system mounting length can be calculated. 
This mount length is the distance between the bottom of the propellant tank and the gimbal attach 
point at the top of the reactor where the turbomachinery, lines, and valves are located. If a value 
other than 0.0 is input for the mount length, XMOUNT, that value will be used in finding the total 
engine length. If XMOUNT = 0.0, it is calculated by the program as: 

XMOUNT = 0.5*Reactor length + Total pump length 


where. 

Total axial pump length = NstagefttWDfDia)^ ^ + «L/D)*Dia )lrtjai , + 

<(l7D)*Dia )nirt(jne 

Total centrifical pump length = ((UD)*Dia)^„ ^ + ((UDl'Dia)^ + 


VD = length to diameter ratio (per stage for axial pumps) 

Dia = diameter 

An expander engine cycle is considered balanced when the ratio of required pump 
horsepower to delivered turbine horsepower is approximately equal to 1.0. If the cycle is not 
balanced, a new value for turbine pressure ratio is calculated and the entire design process is 
repeated. For the gas generator and bleed cycles, the turbine mass flowrate is calculated based on 
the horsepower required by the pump and boost pump/inducer, and a balance is achieved in this 
manner rather than through pressure ratio adjustment. 

Some important, yet easily overlooked inputs are fluid specific heat ratio, GAMGPB, and 
heat capacity, CPGGPB. Despite the misleading variable names that seem to refer only to the GG 
cycle, these values are used for all cycles in various locations, such as the turbine enthalpy 
calculation. In general, the code calculates the heat capacity at each point where it is needed based 
on conditions at that point, but CPGGPB is used often enough to require a reasonable value be 
input, and it can be the factor that causes an expander cycle to either balance or fail. Default values 
for these variables are for a mixture of hydrogen and oxygen at a ratio of 0.75. However, for code 
operation using an expander or bleed cycle, these values should be set to values appropriate for 
hydrogen at similar condidons as those experienced at the turbine inlet during the cycle of interest. 
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Some sample values for hydrogen at 140()°R and 1000 psi are GAMGPB=1 .46 and 
CPGGPB=3.51; these values were obtained from the new hydrogen properties data (see Ref. 2-7). 

An important input for expander cycle TPA design is the turbine bypass ratio, BYPTUR; it 
is the ratio of reflector outflow that goes directly to the core divided by the total reflector outflow. 
The tie tube flow goes directly to the turbine and is therefore not affected by this bypass. As the 
bypass ratio acts only on the reflector flow, the user must be careful when determining this value. 
For example, if an overall turbine bypass of 50% is desired and the nozzle flow fraction is 0.70 
(30% of flow goes to tie tubes, 70% to nozzle), the turbine bypass ratio BYPTUR is calculated and 
input as 0.5/0 .7 = 0.7 1 . Setting BYPTUR equal to 1 .0 will cause all reflector flow to be dumped 
directly into the core so that the turbines are driven by tie tube flow only. 

The gas generator cycle requires input of the GG mixture ratio, OFGGPB, the ratio of 
specific heats, GAMGPB, the specific heat, CPGGPB, and the molecular weight, WMGGPB. 
The default values for these variables are for LOX/H 2 at approximately 1400 psia. The ratio of 
specific heats, specific heat, and molecular weight were determined by a run of the ODE module of 
the TDK computer code using the desired pressure and mixture ratio. The user also inputs the 
turbine outlet pressure, PTURBO, and the pressure ratio across the gas generator/pre-burner, 
PBPFR, PBPRO. For the bleed cycle, the user selects the analysis method and then inputs either 
the turbine inlet temperature, TURBTIN, or the bleed mass flow fractions, FRACHB, FRACCB, 
along with turbine outlet pressure, PTURBO. 

The multiple propellant feed leg TPA option (JCNFIG=5) was added to ELES to allow for 
the redundancy usually desired with NTP engines. Typically, two feed legs will be desired, with 
one half of the total flow running through each pump and turbine during normal operation, as can 
be seen in the cycle schematics in Figure 2-3. This option is normally used with the pump-out 
(double run) option as described in Section 2.3.1. If three feed legs are desired, the initial pump- 
out design run is made assuming one pump is out and two are operational. These two remaining 
pumps are typically designed to handle the full thrust level (FFRAC = 1.0). When multiple feed 
legs are used, the TPA output lists the weight for each pump and turbine in their corresponding 
output sections, while the final engine system TPA summary section lists weights for the total 
turbopump feed system. 

Another new code option is the evaluation of a user-defined TPA, which is described in 
detail in Section 2.3.3. This option allows evaluation of off-design pump and turbine 
performance. It is used automatically with the double run option in which turbomachinery is 
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designed at a pump-out thrust level and then multiple pumps and turbines possessing the 
previously determined characteristics are evaluated at full thrust level. The flag to initiate the user- 

defined TPA design option is ISTSET = 1, and INPTPA = 1 indicates that TPA-related weights 
will be input. 

2.2.11 Weight Multipliers 


Due to the wide range of possible design strategies available for most engine components, 

weight multipliers are provided for all major components. These multipliers are useful when trying 

to match existing designs or design methods. They are also used to account for excess component 

weight not specifically calculated in the code; for example, the standard tank weight multiplier is 

1.7 to allow for the extra material required for weld lands and fittings, see Ref. 1-2. Some of these 

weight multipliers have been discussed in detail elsewhere in this report; all will be summarized 
here. 


The weight multipliers are listed in the worksheet, see Appendix A, along with their default 
values. All tank-related multipliers are set to 1.0 as NESS will primarily be used for engine 
design; the user must input any desired value other than this default. The total nozzle and haidware 
multiplier, CXWENG, is set to 1.0 as it is more likely that the multipliers for individual 
components will be used to account for extra weight rather than adjusting the entire engine weight. 
The valve multiplier, CXVALV, is set to 2.8 to account for dual valves (for redundancy) and a 
factor of 1.4 to include some extra valve weights (other than the main valve) not explicitly 
calculated in NESS. The convergent nozzle multiplier, CXWCHM, is set to 1.0. CXWNZE is the 
nozzle extension multiplier and is used on all portions of the nozzle extension (tubes + radiation- 
cooled portion when used); its value of 1.1 allows for flanges and fittings. 


Hot gas ducting weight is adjusted with CXWDUC that is set to a value of 3.5 to account 
for the weight of flanges, bolts, bellows, bosses, insulation, etc. The gimbal system (excluding 
the power supply) is multiplied by a factor of 1 .4 as set by the variable CXWGIM. The thrust 
mount multiplier CXWTHM is set to 0.9 to allow for technology advances not included in the 
NERVA-era weight correlation between thrust structure and reactor power. The gas generator 
injector weight is multiplied by 1.4 as input by CXWIGG. The turbine weight is multiplied by a 
factor of 1.3 using CXWTPA, and all pump weights are multiplied by CXWPMP, a factor of 
either 1.3 (centrifugal pumps) as was deemed necessary after comparison with other engine 
designs. The multiplier for axial pump weights depends on the thrust level, with a value of 4.93 
recommended for thrust levels below 50,000 lbf, a value of 5.75 used over the range 50,000- 
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100,000 lbf, and a value of 6.0 used for thrust greater than 100,000 Ibf. A comparison with the 
few existing design weights, and the multiplier used to achieve these weights is displayed in 
Table 2-3. 


Table 2-3. Axial Pump Weight Comparison and Multipliers 



Single Lea 1 

Dual 

J £1 J 

1 


■TQFfEZS 

■v'l 

hudh 


■m 

25,000 

90* 

87.8 

130* 

141.7 

4.93 

50,000 

200* 

172.7 

270* 

274.0 

5.75 

75,000 

— 


400* 

396.9 

5.75 

104,000 

1030** 

1037.5 

— 

— 

6.0 


*ftef. 2-9 “ 
**Ref. 2-10 


Comparison with existing designs gives an ignition system multiplier CXWIGN with a 
value of 1.3. Engine bay lines are multiplied by 2.5 to allow for flanges, bolts, bellows, etc. The 
TPA components, valves, and engine bay lines are all multiplied automatically by the number of 
propellant feed legs, NTPA, when appropriate. 

The support hardware multipliers, CXWPNEU, CXWINST, and CXWTNKAS, are 
discussed in the support hardware section of this report, and reflect the technology advances made 
since the correlations used to calculate the component weights were developed. 

2.3 Additional Features 

A n um ber of features have been added to the original ELES to more accurately model a 
nuclear thermal propulsion system. 

2.3.1 Pump-Out Option 

A typical nuclear propulsion system will include multiple propellant feed legs for 
redundancy. Each feed leg will be designed to a desired pump-out thrust level that is less than the 
nominal operating value. To accurately model this feature, a computer run would have to be made 
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al this reduced thrust level to design/sire a single pump and turbine for these conditions, and then 
these values would be used for a second run at full thrust level with multiple pumps to determine 
nominal operating conditions. To simplify this process, a pump-out (double-run) option is 
available for ail engine cycles. The first pass through the code designs a single shaft turbopump 
that operates at a reduced thrust level and corresponding reduced chamber pressure (pump-out 
condmons) specified by the user. The second pass automatically assigns the pump and turbine 
parameters calculated by the first run to be inputs for the user-defined TPA option. The valve and 
engtne bay line weights from die first run are also retained to be output with die total engine 
summary. The second pass will design a system using an input number of identical pmpellant feed 
legs, each with characteristics as calculated in the first pass. 


If the number of feed legs is greater than two. the first pass through the code, the pump-out 
run. will be made assuming NTPA-1 feed legs. For example, if the desired number of feed legs is 
3, the first run will be made assuming a single pump failure and will perform the cycle analysis 
with 2 pumps. The individual pump, turbine, valve, and line weights will be retained from the 
first run and later multiplied by the total number of feed legs as usual. 


To utilize this option, the input file must contain IDBLRUN = 1 and a corresponding thrust 
Wei fraction FFRAC (default = 0.8, or 80% thrust level), based on a dual propellant feed system. 
The user must set the pump configuration flag to the single shaft option, or JCNFIG = 2; the code 
automatically sets JCNFIG = 5 and assigns the pump and turbine parameters calculated in the first 
pass to the appropriate user-defined TPA variables for the second pass. In the input file, the user 
specifies the number of identical feed legs to be used for the second pass as NTPA 


Upon completion of a double run, the user must examine the turbopump output for the off- 
design run to detennine whether the design is feasible for operation in both thrust regimes An 
opuon has been added to NESS to be used with the double run option that will perform this 
ev uation automatically, and if the design fails either of the tests used, the thrust fraction FFRAC 
set for the in, ual pump design run will be reduced by 5% and the entire process repeated until either 
an acceptable design ,s achieved or the thrust fraction becomes less than the fraction of flow 

IrRA^T h ^ UmP ^ fU " thrUS * ,eVd (i CM FFRAC < 1/NTPA) - To sclect «“* °P tion > 

ITRATE=1. The tests used to determine adequate off-design performance include a check for axial 
pump specific speed above 800. Also, the axial pumps cannot be throttled below 60% and 
therefore a test ,s made to determine whether the axial pump volumetric flow is at least 60% of the 
metnc flow handled by the pump in the initial low-thrust design run. This same test is 
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performed for the centrifugal pumps, with the throttling limit set to 40% instead of the 60% limit 
used fra* the axial pumps. 


For all runs made using the pump-out option, whether iterative or not, the speed (RPMs) of 
the pumps calculated during off-design (full thrust) operation is compared with the calculated blade 
root stress speed limit. This speed limit is based on turbine size and blade material properties. If 
the calculated pump RPMs are more than 3% higher than the speed limit, a warning is pnnted out 
in the warning section of the output indicating that pump RPMs are too high and pump design is 
nominal. This speed limit problem can often be overcome by using a higher pump-out design 

thrust fraction (FFRAC). 

2.3.2 User-Defined Engine Burn Time 

An option has been added which allows the user to input the engine bum time rather than 
have the code calculate the bum time based on flowrates and input amount of propellant. This 
option is useful when the amount of propellant to be used is unknown or the tankage design is not 
important. This bum time is used mainly to size the gimbal power supply, whose weight is time- 
dependent. To use this option, set the flag IUSRBRN equal to 1 and then input bum time in 

seconds as TUSRBRN. 

2.3.3 User- Defined Turbomachinery 

The user-defined turbomachinery option of NESS allows evaluation of pump and turbine 
performance at off-design operating characteristics and with a variety of propellants. The 
parameters input to define the TPA for off-design evaluation are detailed in the worksheets 
following, and include number of stages for all pumps and turbines, pump and turbine diameters, 
turbine annulus area, turbine admission fraction, and various gas generator/mixer parameters. 


NESS calculates pump head rise and volumetric flowrate, and turbine horsepower, mass 
flowrate, and pressure ratio based on cycle balance requirements. For the centrifugal pumps, these 
values are used to calculate the pump rpm as a function of input pump diameter. To perform this 
calculation, a correlation had to be developed for pump head coefficient as a function of specific 
speed (sta ndar d cases interpolate this coefficient from a data table), and is of the form. 


HC = const * SS X 
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where 

HC = head coefficient 
SS = pump specific speed 


For example, the main pump correlation is: 


HC = 3.7852 * SS-0-28786 


This correlation is different for main pumps and boost pumps. The specific speed is a funedon of 
pump rpm, head rise, and volumetric flowrate, as is shown below: 

SS = RPM * SQRT(volumetric flowrate)/(pump head rise°-75) 

The pump diameter is calculated as: 


Dia (720/pi*RPM) * SQRT(32.2*pump head rise/head coefficient) 

Substituting the head coefficient and specific speed equations into the equation for pump diameter 
and rearranging gives an equation for pump rpms as a function of input pump diameter only Once 

the rpms are known, the specific speed, efficiency, and horsepower are easily found from the 
standard ELES equations. 


The axial pump user-defined TPA method is slightly different from that used for the 
centrifugal pumps. Using the equations listed above for specific speed (SS) and pump diameter 

along with the best-fit equation determined for pump head coefficient as a function of specific 
speed, the following equation is found: 


52525*32.2*QFL/(Dia 2 *sqrt(HFL))=SS 2 (0.88237-2.3145E-4*SS+ 

2.3161E-8*SS 2 -7.7028E-13*SS 3 ) 

where: 

QFL = volumetric flowrate (gpm) 

HFL = pump pressure-head rise (ft) 

Dia = pump diameter (in.) 


This equation is a function of specific speed only, and is of the fotm of a fifth-order polynonual 
that can now be solved by the iteradve secant method. Once the specific speed is known, the pump 
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rpms can be found and the rest of the calculation proceeds the same as for the centrifugal pumps. 
The inducer is again modeled as a boost pump with its speed (rpm) and diameter fixed to that of the 
main pump. The inducer head rise is determined by the input fraction BPFRFL, and the specific 
speed, efficiency, and horsepower are now easily calculated by the standard equations. 


The user-defined TPA option of NESS calculates the required turbine mass flowrate and 
horsepower and then evaluates the user input turbine to see how well it performs in meeting these 
requirements. The first step is to calculate the isentropic spouting velocity (Co) based on the 
number of turbine stages. Now calculate the ratio of turbine blade tangential velocity to Co based 
on input turbine diameter (U/Co) and check whether this ratio is within the accepted range of 0.2 - 
0.6; if not, print a warning. Next, calculate the turbine inlet mach number and check whether it is 
below the accepted maximum value of 1.7; issue a warning if not. Finally, calculate turbine 
specific speed, efficiency, and horsepower provided. Compare the horsepower provided with the 
horsepower required and if not within 3%, calculate a new turbine pressure ratio and repeat the 

entire process. 

To use this option, first set the variables ISTSET = 1 and INPTPA=1 to indicate that the 
TPA is user-defined and the TPA-related weights will be input. The number of pump stages are 
input with PDIAFL and PDIAOX. Turbine stages are input with either TSTGES for a single shaft 
turbine, or TSTAGF and TSTAGO for fuel and ox turbines (can be used only for GG cycles). 
Diameters are input in inches with PDIAFL and PDIAOX, and either TDIAM or TDIAFL and 
TDIAOX. Boost pump diameters can be input with BPDIAF and BPDIAO. Turbines also need to 
have admission fraction and annulus area input using the variables listed in the worksheet. TPA- 
related weights will not be calculated for the user-defined TPA option and therefore the user may 
input these weights for total TPA, TP AWT, start system, WSTART, ignition system, WIGNIT, 
hot gas manifolding, WHGMF, autogenous heat exchanger, WTHTX, and gas 
generator/prebumer, WGGPB. If not input, the weight summaries will list these weights as zero 
unless a double run is being made, in which case the weights calculated in the first pass are 
retained and printed out in the output summaries. 


The user-defined gas generator cycle requires many more inputs than are required for the 
expander cycle. First set the flag IUSRGG equal to 1 to indicate a user-defined GG and input all 
pump and turbine parameters as described above. In order to insure that the GG and turbine are 
modeled correctly, the turbine inlet and outlet pressures, PUSRTI and PTURBO, respectively, 
must be set to the values calculated/input for the NESS-calculated case. For example, if a NESS- 
calculated GG cycle using LQ2/H2 is designed at 80% thrust level and is next to be evaluated at 
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50% thrust level, the turbine inlet and outlet pressures calculated by NESS in the first run must be 
used as inputs for the user-defined run. The turbine inlet temperature, TUSRGG, should be set to 
the actual value found for the propellant combination at given mixture ratio and pressure; normally 
this temperature will simply be the same as that found in the 80% run. If a different propellant is to 
be evaluated or the GG is being input based on an existing design (not NESS-generated), this 
temperature can be found most easily by an initial NESS run where the user-defined option is not 
used and the GG is at conditions similar to those to be used for the actual user-defined run. The 
turbine flowrate, although listed as an input, is actually calculated by NESS as the correct amount 
of fluid flow required for the given operating conditions. The GG bleed flowrate, Isp, and 
efficiency can be set to any reasonable values. The drive fluid parameters must be input to any 
value other than zero. 

Inputs required for the user-defined bleed cycle include IUSRGG=1 and a matched turbine 
outlet pressure PTURBO and inlet temperature TUSRGG. The standard pump and turbine 
parameters, such as diameters and number of stages, are retained or input as usual. 

2.3.4 Weight Margin 

The user may now input a fraction of the total non-nuclear weight to be added in as a 
margin weight. Inside the code, non-nuclear weight is the sum of nozzle weight, total TPA 
weight, lines, valves, thrust mount, support hardware, and total gimbal system. The percent 
(fraction) of this weight to be used as margin is input with FMARG, whose default is 0.02 (2% 
margin). In the output summary, the "non-nuclear weight" includes the weight margin. 

2.3.5 Bleed Cycle Component Models 

The bleed cycle requires a number of extra lines and valves which differ from those 
required by the other engine cycles. The hot bleed flow is tapped off the reactor chamber flow at 
the exit region. NESS assumes that a single line is used for all hot bleed flow regardless of the 
number of TPA propellant feed legs. If required, the hot bleed flow is split into the necessary 
number of feed lines prior to entering the mixers. Cold bleed flow is tapped off from the pump 
outlet flow, and one line is needed for each feed leg. Each cold bleed line includes a valve that 
steps the pressure down as needed to meet the required mixer inlet pressure. One mixer is used for 
each propellant feed leg, and consists of a hot bleed line wrapped by a cold flow line. The flows 
are merged at the mixer outlet and then sent to the turbine. Each turbine has its own inlet line 
containing a throttling valve for flow regulation. Pressure drops across all lines and the turbine 
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throttling valve are input as fractions of the line/valve inlet pressure using the variables CPLINH, 
CPLINC, CPLINT, CPVLVT. The cold bleed valve pressure drop is calculated automatically 
based on the mixer pressure requirement. The hydrogen velocity was assumed to be 200 ft/sec, a 
typical value used elsewhere in NESS for line calculations. For pump-out (double-run) cases, the 
bleed components are designed on the second pass (the full thrust run). 

Table 2-4. Bleed Line Component Design Characteristics 



Equivalent Length 

Flange 

Bellows 

Elbow 

0.2*Line diameter 
1.0* Line diameter 
5. 5* Line diameter 


Table 2-5. Bleed Cycle Line Characteristics 



Number 

Number 


Comoonent 

of items 

Flanges 

Bellows 

Elbows 

Line Lenath 

Cold Need line 

NTPA 

2 

2 

2 

1.5*Reactor Pressure Vessel 
Diameter 

Hot bleed line 

i 

2 

1 

2 

1 .5*Reactor Pressure Vessel 
Length 

Turbine inlet line 

NTPA 

2 

2 

3 

0.7*Reactor Pressure Vessel 
Diameter 

Mixer 

NTPA 

2 

1 

0 

3.5*Turbine Line Diameter 


Each flow line consists of the line itself plus a number of flanges, bellows, and elbows. 
Using a method described in the TRW report, see Ref. 1-1, each line length was calculated as a 
series of equivalent lengths, most of which are functions of the line diameter. It was assumed that 
each flange would add a length that was 20% of the line diameter, a bellows adds 100% of the line 
diameter in length, and an elbow has a length 5.5 times the line diameter (sec Figure 2-4). The 
lengths of the lines themselves are functions of the reactor diameter and length. A summary of 
bleed line component design characteristics is shown in Figure 2-5. Total assumed line lengths are 
as follows: 
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Cold bleed line = 1.5‘D^^ + 2 bellows + 2 elbows + 2 flanges 
Hot bleed line = l-5*L reactor + 1 bellows + 2 elbows + 2 flanges 
Turbine inlet line = 0.7*D reac(or + 2 bellows + 3 elbows + 2 flanges 

Line diameters are calculated as a function of mass flowrate, and fluid density and velocity. 
Line thickness is a function of pressure, line diameter, and line material strength, see Ref. 2-1. 
Once the line length, diameter, and thickness are known, the volume of material can be found and 
finally the weight of each line is determined. Hot bleed line weight assumes all hot bleed flow 
travels through a single line, while each cold bleed and turbine inlet line are sized assuming 
( l/NTPA)*flowrate. 

The cold bleed and turbine throttling valves are sized using the standard NESS/ELES valve 
weight procedure. Valve weights are a function of valve material density, mass flowrate, pressure 
drop across the valve, and fluid density in the form: 

Valve Weight = 1.476*rho v j v *mdot/sqrt(rho H2 *delta-P) 

where: 


rho vlv 

valve material density (lb/in 3 ) 

mdot = 

fluid mass flowrate (lb/sec) 

rh°H2 = 

fluid density (IbAn 3 ) 

delta-P = 

pressure drop across the valve (psi) 


The typical valve multiplying factor, CXVALV, of 2.8 (=2*1.4) is applied to the bleed line valves 
along with all other valves in the cycle. The factor of 2 doubles the number of required valves to 
provide redundancy. For example, a typical bleed cycle has two feed legs and therefore requires 
one cold bleed and turbine valve for each leg; application of the valve multiplier will instead allow 
for two valves of each sort for each feed leg to satisfy the usual redundancy requirements. The 1.4 
factor accounts for valve weight not specifically calculated by NESS, see Ref. 2-1. 

The mixer design is shown in Figure 2-11. It consists of diverging and cylindrical cold 
flow portions plus a hot flow line. The overall mixer length is assumed to be: 
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Figure 2-11. Bleed Cycle Mixer Design 

Mixer Length = 3.5*D turbine , nlet + 1 flange (cold side)+l flange (mixer outlet) 

+ 1 bellows (mixer outlet) 

where: 

D tuibine inlet = Turbine inlet line diameter 

The hot flow line length within the mixer is assumed to be: 

Hot Line Length = 2.5‘D^^ inlet + 1 flange + 1 bellows 

The length of the diverging section of the cold flow portion of the mixer is set to: 

Divergent Length= 2*0^^ mUtl + 1 flange 

The volume of material needed for each portion of the mixer is now calculated individually so that 
the weight for each portion, and finally the entire mixer weight, is determined. 
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2.4 Code Setup and Execution 


NESS is written in FORTRAN 77 and currently resides on a VAX mainframe computer 
system. It has recently been modified for use on a personal computer (PC) as well. Allowing a 
much larger range of potential users. The major NESS modification required for transport to a PC 
was a reformatting of the input file from a VAX-specific namelist file to an unformatted read file 

that contains all possible input variables along with a brief description of each variable's function, 
see Appendix A. 

The entire NESS code is made up of four parts: the source code, the executable, the library 
of subroutine object files, and a library of propellant performance data. The NESS executable 
takes up approximately 4700 blocks of storage space on the VAX and about 1 .6 MB on a PC. The 
source code for NESS is made up of approximately 220 subroutines that have been separated into 
individual files for easier editing. The subroutines take up approximately 4800 blocks of storage 
space on the VAX, and 1.9 MB on a PC. The object library ELES_LIB.OLB takes up about 
6,000 block of storage space, or about 2 MB on the PC. The propellant performance library is 
included with the code, but may not be needed as all hydrogen performance data has been entered 
elsewhere in the code; this data uses 72 blocks of storage. If storage space is a problem, the 

executable alone could be loaded onto the computer while the rest of the code is left on disk to be 
loaded as needed. 


The standard NESS operation requires creation of a structured series of directories 
whenever the code is loaded onto a new computer system. The executable and propellant data file 
must be put into a directory called [account name.ELESJ. The input files reside in the directory 
[account name.ELES.INPUT], and the output appears in the directory [account 
name.ELES. OUTPUT]. The source code and object file library are loaded into [account 
name.TEMP.CURRENT], If the code will be run in debug mode, a directory [account 
name.TEST] must be set up that includes the propellant library file PROPLIB.DAT and an input 
file with the name ELES.INP. This directory structure is especially suited for VAX operation. 

To simplify the code structure, PC users may wish to edit the governing .COM (or batch) 
files - ELES_SETUP.COM and RUN.COM - to allow placement of all code parts (executable, 
source, input, output, etc.) into a single directory. The code itself can be easily edited to operate 
off a single input file (always the same name), if desired. NESS has been successfully tested on a 
PC using Lahcy FORTRAN 77/EM-32 with the entire code contained in a single directory. The 
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code was edited to always open a certain input file and can be run by simply typing the word 
NESS (or whatever the code file is called); .COM files are not required at all. 

For standard (VAX) code operation, a number of *.COM files are necessary/useful for 
code execution. The file ELES_SETUP.COM must be run at some point before the code is run to 
insure proper directory and file initialization; this is most simply achieved by adding this file to the 
LOGIN.COM file and having it execute automatically with each login. In the [..CURRENT] 
directory, the file FL.COM is used to compile an individual subroutine and add/replace it in the 
object library; it is used as "@FL filename ". FALL.COM will recompile all subroutines and 
replace their previous versions in the object library. To link the governing routine with the object 
library, type "@LD" to execute LD.COM and a new NESS executable will be created. 

If the code has been edited to always open the same input file, all input files must have the 
name assigned by the programmer (NESS.INP, for example), otherwise, the filename must have 
the extension .inp and must contain 10 characters or less, excluding the extension. To run the code 
(standard operation, unmodified code) type "MODEC filename" without the filename extension of 
.inp; for example, typing "MODEC NTPREGEN" will run NESS with the input file 
NTPREGEN.INP and place the output in a file called NTPREGEN.OUT in the output directory. 
A file called NTPREGEN_ELES .OUT is also created in the output directory that is essentially a 
printout of the input file. If the computer has a debug mode, enter the [account name.TEST] 
directory and type "RUN ELES: MODEC" and the code will execute using the input file 
ELES.INP. For PC operation, the programmer/user may compile and link all source code into an 
executable, called NESS for example, and edit the code to open a single input file, such as 
NESS.INP. If this is done, the user need only type "NESS" to ran the code, and the output will 
appear in the same directory as the code with a name assigned by the programmer. 
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3.0 REACTOR SYSTEM 


scaling relations. 


internal shield, modeling assumptions, and 
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3.1 Reactor System Descript 

An engineering description of the reactor's major subassemblies for both the ENABLER I 
and II reactor systems am given in the following secdons. ENABLER I 

3.1.1 Reactor Assembly 

formailtllTr ^ - - actuation system 

or reacttvtty contiol devtces wtth associated insmtmeniation and controls am shown in figure 3 - 1 

reactor consists of fuel elements, support elements, a com periphery, suppori plants aliens 

an mtema s teld, a reflector assembly, and control drum drive assemblies. Reflector coolant is 

zZT f :z 

routiuHr n PUmP ' Addid0nal ,urbop “ m P flow ™y also be obtained by 

r rr ,0r C °° lam ' XhaUS ' » * "AW ^ ttubine exhaust gas flows through , e 

^“rrr- p,ena i -— * - ^ p - - i szz 

“ A 8 reaC ‘° r ™ S *“ is to ‘ Bd "y ft' reactor assembly to operatic 

temperatures and exhausted out the nozzle. ^ g 

3.1.2 Fuel and Support Elements 


The fuel elements in Figure 3-2 for the FNARi pd „i„~ 
flinr . t . , 6 ‘"“e tiNABLER class reactor serve the combined 

excha er su^ r^e 7” ^ ^ * hydroe '" “ d - » - 

elen« T ^T , V n ? y “ Pr ° V ‘ d ' d ‘ h " >U8h flSSi °" ° f ^ contained in the fuel 
ent. Table 3- 1 lists the characteristics of the three fuel materials defined in the NESS code 
Multiple coolant channels coated u/ith VrC' tr ■ . c code, 

through the elements Th ( for graphite and composite) form flow passages 

also £7 n ' X, ' n ° r ° f th ' h ' XagOM ' fu " *— <«cep, carbide, are 

Thefuel I a s coanng protects the carbon from reaction with the hydrogen propellant 
The fuel element dimensions were established by the NERVAnn™, • propetlam. 

hexagonal with 19 holes 0 100 to 0 1 10 inch 7 " ' a nom,nal 0 75 

coating thicknesses « se,„“ n ^ ^ df'", T" ^ - 

^ ^ default values for user variables The rr\A 

pernnis ihe user to specify scaled fuel see Refs 3 I and ti a , SS COde 

allowable fuel power density. <f 3-2. tbai allows an increase in the 
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Figure 3-1. ENABLER Class (NERVA Type) Nuclear Thermal Rocket Engine Reactor 
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Table 3-1. Fuel and Support Element Parameters 


Fuel Element 
Composition 

Graphite 

Composite 

Carbide 

Temperature Range (°K) 

2200-2500 

2500-2900 

2900-3300 

Fuel 

Coaled Particle 

UC ZrC 

Solid Solution and Carbon 

(UZUC 
Solid Solution 

Coaling 

ZiC 

ZrC 

— 

Unfueled Support Element 
Composition 

Graphite 

ZiC -Graphite Composite 

ZiC 

Unfueled Element Coating 

ZiC 

ZrC 

1 


The reactor core is supported in the longitudinal direction by tie-tubes anchored in the core 
support plate above (inlet end) the cote. The de-tubes run the full length of the reactor core and 
connect to fuel support blocks below (exhaust end) the core. These de-tubes are located inside 
unfueled support elements, which have the same length and external dimensions as the fuel 

elements. Hie support elements have a single, large longitudinal hole with a P^^ nsula “« 
liner. Within the hole is located the tie tube assembly, which may incorporate ZrHj moderator 

required. The support element composition is given in Table 3-1. 


The reactor core is sized based on an average fuel element power of 1.2 MW per element 
and one support element per six fuel elements in Table 3-2 at thrust levels greater than 50,000 
pounds. The 1.2 MW per fuel element was demonstrated in the Pewee reactor (402 fuel elements 
with a power level of 503 MW) and was the design level for the Phoebus-2A reactor (4068 fuel 
elements with a 5000 MW design power level). For the smaller reactors, sufficient reactivity is 
obtained by increasing the relative number of support elements to fuel elements (Table 3-2) which 
increases the amount of rirconium hydride moderator to the desired level. Also to keep a 
reasonable core length to diameter ratio (<2) for the smaller reactois (15000-25000 Ibf. thrust) the 
element length was set at 35 inches. At the 25000 thrust level (Pewee sire core volume) the 
relative power density of the fuel element is the same as the larger reactors (1.2 MW/52 inch). 
However, at the lowest thrust level (15,000 Ibf.) the fuel element power density had to be reduced 
in order to obtain a core large enough for criticality. Fuel volume, xireonium hydride loaduig, and 
reflector thickness all act to increase core reactivity. Neuoonic analysis to determine the exact 
axMtation of these parameters that achieves criticality is not part of the NESS code at this time. 
The NESS code does provide a warning message if the selected combination of parameters ts 

questionable. 
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Table 3-2. Reactor Parameters as a Function of Thrust Level 


Thrust (lbf) 

Reactor Power Range 

Fuel and Support Element Length (inch) 

Pressure Vessel Length (inch) 

Fuel Element Power (MW) 

Relative Fuel Element Power Density 

Ratio of Fuel Elements (N) to Support Elements 


15,000 

25,000 

>50,000 

275-400 

460-670 

920-6700 

35 

35 

52 

82.6 

84 

101.6 

0.629 

0.808 

1.20 

0.778 

1.0 

1.0 

2:1 

3:1 

6:1 


3.1.3 Radiation Shield 

A radiation shield internal to the pressure vessel is used to reduce the gamma and neutron 
flux levels in the engine components forward of the reactor. This internal shield limits radiation 
leakage through a plane 63 inches forward of the core center, perpendicular to the engine axis, to 
the levels given in Table 3-3. The shield is located immediately upstream of the core support plate, 
see Figure 3-1. The reactor internal shields for the thrust levels over 50,000 lbf. have about 12.5 
inches of Borated Aluminum Titanium Hydride (BATH) and about 1.3 inches of lead. At the 

lower thrust levels the thickness of the BATH and lead is slightly reduced due to lower core power 
density. 


Table 3-3. Radiation Leakage Limits at a Plane 63 Inches Forward of the Core Center 


Type of Radiation 

Radiation Leakage Limits 
Within Pressure Vessel 
Outside Radius 

Gamma Carbon KERMA Rate 

1.8 x 10 7 Rad(c)/hr 

Fast Neutron Flux 

2.0 x 10 12 n/cm 2- sec 

Intermediate Neutron Flux 

3.0 x 10 12 n/cm 2 -sec, 
0.4 eVsEnS 1.0 MeV 

Thermal Neutron Flux 

6.0 x 10 n n/cm 2 -sec 
En < 0.4 ev 
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3.1.4 Reactor Propellant/Coolant Circuits 

In an NTP system, a nuclear reactor supplies the energy to heat the propellant flowing 
through the engine. The hot propellant flows into a nozzle that functions in the same manner as a 
chemical engine. The reactor in an ENABLER reactor-based NTP engine system generally has 
three propellant (coolant) circuits as shown in Figure 3-3. The primary circuit is through the 
central shield and core into the chamber. This circuit provides more than 90% of the heat to the 
propellant. All the components surrounding the core require cooling due to the radiation induced 
heating and heat transfer from the primary stream. The propellant cooling of the ex-core 
components is divided into two additional circuits: the tie tube (core support) circuit and the 
peripheral component circuit that includes the core reflector and extension shield. These circuits 
along with the nozzle regenerative cooling circuit provide the first pass through the reactor system 
few the propellant, which acts as component coolant The heat supplied by these secondary circuits 
provides the energy to power the turbopump. After passing through the turbine, all the propellant 
passes through the primary core circuit and into the nozzle to provide the engine thrust. 


COMPONENT BLOCK DIAGRAM 



HEAT GENERATION 

Coro 

-1,500 MW 

Tie Tubes 

3-7% 

Reflector 

1-2% 

Central Shield 

-0.2% 

Ext. Shield 

-0.03% 


Figure 3-3. Propellant Flow Circuits Through the Reactor 
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The NESS code allows the user to choose one of two propellant circuit connection options. 
The first option routes the outlet of both the tie tube circuit and the peripheral component circuit to 
the turbopump. This arrangement was selected for the NERVA flight engine, the R-l, and it 
provides sufficient energy to the turbopump to allow operation of the engine at high chamber 
pressures (approximately 2000 psia) for small engine systems (<40,000 lbf thrust). The second 
option routes only the tie tube circuit outlet to the turbopump, while the outlet of the peripheral 
component circuit is routed directly to the core inlet. This configuration saves weight by 
eliminating a massive flow baffle at the top of the core and by reducing the pressure vessel design 
pressure, thus decreasing its thickness. With this configuration, the energy available to drive the 
turbopump is reduced and therefore engine operation at high chamber pressures ( > 1000 psia ) 
may not be possible in engines with a 6:1 fuel to support ratio. The second configuration is 
generally preferable if a cycle balance can be achieved. 

3.2 Baseline Reactor Design 

The Rover/NERVA database provides numerous reference designs for reactors and engines 
in the size range of 15,000 lbf to greater than 250, 000 lbf thrust range. The engine modeled in the 
NESS program is the ENABLER reactor class of NTP engine systems, which is discussed in Ref. 
1-4, that is derived from the nuclear rocket technology developed in the Rover/NERVA programs. 
The ENABLER designs incorporates NERVA type fuel elements which are 0.75 inch (19 mm) 
hexagonal extrusions of graphite based fuel with a 19 coolant channel array within the element. 
The code allows the user to select from one of the three fuel materials developed during the 
Rover/NERVA program: Graphitic, Composite, or Carbide. The ENABLER engine is generally 
specified with fuel elements fabricated from the (U.Zr)C-Graphite composite material developed 
late in the Rover/NERVA program, which exhibits improved corrosion resistance and allows 
higher operating temperatures and power densities, see Refs. 3-3 and 3-4. Zirconium-hydride 
moderator is placed in the core support elements (demonstrated in the Pewee reactor) to increase 
the neutronic reactivity and thereby decrease the required uranium fuel loading. 

Detailed data is available on the breakdown of actual reactor system component masses. In 
the NESS model the core size is based on the number of fuel elements needed to meet the required 
power level. The design of the reactor peripheral regions follows the R- 1 engine design shown in 
Figure 3-4, but the peripheral components are sized according to the core dimensions. For the R-l 
reactor shown in Figure 3-4, the nominal core dimensions are 38 inch (96 cm) diameter by 52 inch 
(132 cm ) long. The components surrounding the core are sized to satisfy structural and neutronic 
requirements. The major components are the core barrel, reflector, pressure vessel, core support 
plate, flow baffles, and top shields. 


3-7 




Figure 3-4. Layout Drawing of the R-l Reactor 


3.3 Reactor Core Design and Thermal-Hydraulic Model 

The required core power level is determined from the specified engine flow and chamber 
temperature. The core power level and the average allowable heat generation of a fuel element 
determines the total number of fuel elements and support elements in the core. Based on the core 
peaking factor, a single channel analysis is performed to calculate the thermal and pressure profile 
for the peak channel of the peak element in the core. The calculation uses finite increments along 
the channel length beginning at the core exit where the chamber conditions are specified. The 
governing equations are given below. 

The convective heat transfer between the fluid and channel wall is defined by: 
q = h c A $ ( T w - T r ) 

where T w is the channel wall temperature and T r is the coolant gas stagnation recovery temperature. 
For small Mach numbers (« 1 .0) the difference between the recovery temperature (T r ) and the 
fluid free stream bulk temperature (T^) is not significant, so that the equation may be written as: 

q = h c A,(T w -T b ) 

The heat transfer (q) must match the heat generation in the fuel material. The heat 
generation in the fuel is determined by the fuel loading, fuel volume, and neutron fluence. For the 
purposes of the thermal hydraulic calculations it is sufficient to specify a power profile and the total 
power produced by the element. The NESS code uses a cosine power profile typical of that 
observed in the NERVA reactors: 

P = P n cos( 0.89 lit ( x/L - 0.452 ) ) 

where P n is the normalized element power factor and x/L is the normalized axial location in the core 
measured from the inlet. The peak temperature in the fuel (T f ) is determined from the following 

correlations for a heat generating solid with a hexagon array of coolant channels of diameter D and 
pitch S: 


e = 


TtP 2 

3.4641 S 2 
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where k s is the thermal conductivity of the solid. 

The convective heat transfer coefficient, h c , is determined by the McCarthy-Wolf, see Ref. 
3-5, correlation: 


= 0.025-^ Re° 8 
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where the fluid properties are evaluated at the fluid bulk temperature. The entrance effect term 
( 1 + 0.3 ( x/D )‘°- 7 ) is limited to 1.1 for small x. 

As the coolant flows along the channel, it experiences a pressure loss due to wall friction 
and fluid acceleration. The momentum equation for one dimensional flow in finite increment form 
is: 


Pi-P..' = y( v 


1 + 1 


_Vi) + fi fr i(V ' +, + V ‘ ) 


where Pj is the coolant pressure at station i, G n is the mass flow per unit area, v^ is the specific 
volume of the coolant, D h is the hydraulic diameter of the channel, fj is the Fanning friction factor, 
and Dx is the length increment along the channel. The friction factor is obtained from the Taylor, 
see Ref. 3-6, correlation for gaseous flow through a smooth tube: 
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where Re w is a modified surface Reynolds number in which the gas density is evaluated at the fluid 
bulk temperature, but the viscosity is evaluated at the channel wall temperature: 


Re„ = 


fG n D^ 


(V 

V* J 


ItJ 


The evaluation of these equations for the peak channel in the core determines the required 
core pressure drop. 


After the calculation of the core profile and pressure drop, the heat generation rates for the 
core peripheral regions are calculated. Because NESS does not have neutronics analysis 
capabilities, the heat generation in the peripheral regions is defined as a fraction of the total core 
power. After completion of the thermal hydraulics, code control returns to the NESS engine code 
for determination of the cycle balances. 


In addition to the basic thermal-hydraulic parameters of the core, NESS calculates the 
estimated life of the fuel based on the hot end corrosion correlations obtained from the Nuclear 
Furnace 1 and electrical testing, see Ref. 3-3. Fuel life is given by: 


rh = 30.5 exp 


35114 


l l = A f / rh 


\ 

/ 


where T w is the peak wall temperature of the fuel channel at the hot end in degrees Kelvin, is the 
fuel mass loss rate in g/sec per cm of fuel element length, A f is the fuel area, m, jmit is the allowable 
mass loss in g/cm 3 , and t| is the fuel life in seconds. NESS contains the necessary corrections for 
calculating the life of scaled fuel. The fuel life estimate is not valid for the carbide fuel type. 

3.4 Reactor Weight Model 

The reactor mass model divides the reactor system into 53 regions for both types in an R-Z 
model as shown in Figure 3-5 and Table 3-4. Each region contains one, or at most a few, 
components. The masses of all the components and their constituent parts within a region have 
been tallied and converted into a pseudodensity for each region, see Ref. 3-7. The dimensions of 

the regions are based on the core size determined above, with appropriate dimensional dependency 
algorithms. 
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Figure 3-5. R-Z Model of the Regions in the R-l Reactor 
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Table 3-4. Reactor Weight Model Regions 


REGION NUMBER 

REGION DESCRIPTION 

MATERIAL 

1 - 15 

Core 

Fueled Element 
Unfueled Element 
Pyro Sleeve 
A-286 
SS-304 
Hydrogen 

16 

Core Periphery 

Graphitite-G 

Pyrofoil 

ZrC (60% Dense) 
TZM Moly 
Hydrogen 

17 

Lateral Support 

P03 Graphite 
ZTA Graphite 
Pyrofoil 
Hydrogen 

18 

Structure 

P03 Graphite 
AI-6061 
A-286 
Hydrogen 

19 - 21 

Reflector 

P03 Graphite 

Pyrofoil 

Beryllium 

AI-6061 

A-286 

Control Vane 
Hydrogen 

22 

Pressure Vessel Side A 

AI-7039 

Hydrogen 

23 

CHESH 

Pyrographite 

Pyrofoil 

NbC/C Comp. 

W-ThO 

A-286 

SS-304 

SS-316 

Hydrogen 

24 

Nozzle Chamber 

Hydrogen 

25 

Nozzle Barrel 

SS-347 

26 

Aft Reflector Hardware 

AI-6061 

A-286 

SS-440C 

Hydrogen 

27 

Aft Reflector Plenum 

Hydrogen 


TOR29K/14a 
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Table 3-4. Reactor Weight Model Regions (Cont.) 
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Table 3-4. Reactor Weight Model Regions (Cont.) 


REGION NUMBER 



REGION DESCRIPTION 


BATH Peripheral 
Shield IV 


Pressure Vessel Side B 


Lead Central Shield 


Lead Peripheral 
Shield I 


Lead Peripheral 
Shield II 

Lead Peripheral 
Shield III 


Lead Peripheral 
Shield IV 


Peripheral Shield Plate 


Shield Plenum 


Flow Battle II 


MATERIAL 


BATH 

AI-6061 

Hydrogen 


AI-7039 

Hydrogen 


Lead Alloy 

AI-6061 

Hydrogen 


Lead Alloy 

AI-6061 

Hydrogen 


Load Alloy 
Hydrogen 


Lead Alloy 

Hydrogen 

AI-6061 

A-286 

SS-304 

Hydrogen 


Lead Alloy 

AI-6061 

Hydrogen 


A-6061 

A-286 

SS-304 

Hydrogen 


AI-6061 

SS-304 

Hydrogen 


AI-6061 


Central Dome Plenum Hydrogen 


Peripheral Dome Plenum 


AI-6061 

A-286 

SS-304 

Hydrogen 


Pressure Vessel Dome 


Ai-7039 


NERVA Nuclear 
Subsystem 
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The pseudodensity is applied to each region to yield the mass schedule of the reactor for everything 
out to and including the nressure vessel. The weight algorithms automatically delete the flow 
baffles if they are not required based on the choice of flow circuits, see Section 3.1.4. Thrust 
structure, turbopumps, and nozzle masses are not calculated in this module; the NESS code 
determines the balance of engine masses, which is discussed in Section 2.0. 

3.5 Design Variable Options 

User inputs can be divided into three categories: engine parameters, reactor parameters, and 
fuel element parameters. The primary engine parameters are thrust level, chamber temperature, 
chamber pressure, and nozzle expansion ratio. These primary variables are used by the code to 
define the engine specific impulse, propellant flow rate, and required reactor power. The reactor 
parameters include reactor pressure vessel material, power fractions in the peripheral components, 
and tie tube power levels. 

The user supplies the governing parameters for the fuel elements. These include mean fuel 
element power, element dimensions, fuel scaling, and material. The code modules provides for a 
choice from three fuel materials: graphitic (UC 2 beads in graphite), composite ((U.Zr)C-Graphite), 

or carbide ((U,Zr)C). Each fuel type exhibits different properties with regard to mass density, 
power density, and temperature limits. The fuel to support ratio within the core may be set to one 
of three patterns: 2:1, 3:1, or 6:1. The fuel parameters are strictly user defined in that the code 
does not attempt to judge the validity of the inputs. For guidance. Tables 3-1 and 3-2 provide 
information on typical parameters based on the Rover/NERVA technology. The reflector thickness 
and zirconium hydride loading of the tie tubes are also user selected based on the user s estimate of 
criticality requirements. 

3.6 Key Assumptions 

The code assumes that the same basic design will be used at every size level within the 
specified code domain. This provides the basis for calculating the size of the core periphery. 

The code assumes that the user has specified a viable combination of input criteria. For 
example, the code does not verify core criticality and control span. This cannot be accomplished 
until core neutronics is integrated into the code. In particular, engine weights are strongly 
influenced (upward) by criticality considerations for engines with thrust ratings below 50000 lb. 
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Many of the input variables are based on NERVA test experience and should not be altered. 
This includes such things as the limiting fuel element power density and the power distribution in 
the peripheral regions, which arc based on external data sources such as test measurements. 
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4.0 SAMPLE NTP ENGINE SYSTEM DESIGN CASES 


Eight NESS NTP engine design problems are presented in this section. These sample 
design cases demonstrate many of the design capabilities associated with NESS. Key engine 
system design parameters associated with these sample cases are presented in Table 4-1. For each 
sample design case, the NESS VAX mainframe computer input File and engine design output file 
are given (see Tables 4-2 through 4-9). For Sample Case No. 2, initialized NESS program input 
sheets are shown. A clean set of input worksheet forms is given in Appendix A. 
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Table 4-1. Sample Design Case Summary 
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Table 4-2. Sample Case No. 1 
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Input Listing 
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NON-NUCLEAR WEIGHT MARGIN 86.98 
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FLOWRATE SCHEDULE (LB/SEC) FOR STAGE #1 
EXPANDER CYCLE 
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TOTAL ENGINE SYSTEM 

TOTAL ENGINE WEIGHT 12822.2 lb® 5815. 
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GAS PHASE ENCOUNTERED IN RECEN JACKET 

TPA CALCULATIONS TERMINATED BY ACHIEVING DESIRED ACCURACY 
END NOMINAL STAGE DESIGN 




Table 4-4. Sample Case No. 3 
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PRESSURE TANK INSULATION 0.06 
FUEL TANK INSULATION 256.20 
OXIDIZER TANK INSULATION 407.04 
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GROSS IGNITION WEIGHT 34188.44 

GROSS BURNOUT WEIGHT 25131.92 



STAGE DIAMETER 1M.M 

NOZZLE EXIT DIAMETER 148.96 

NUMBER OF NOZZLES 1 

STAGE LENGTH 912.83 
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COLO SIDE FRACTION OF TOTAL BLEED 
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Table 4-5. Sample Case No. 4 


4-118 




Input Listing 


4-119 



*o 


o 

o 

C w 

o — 


*m*-*^m c ~ o 
x T ®» C oo» 
1 *- I I w o- o o 

OUU>«« I 

1>> J*» g«- 
•J JQ U 9 f 
• QQth 


3 0 0 
— 0-0 


®l»OC L 

o • o ^ 
** • o — c 
• — • 
o I o o 
o +* w w ■ 
3 • w 3 « 
O 


:8822i:> 

• o O ** 

t - >(k- 

coo a - 

• f O— ©- c c 5 a o 

# 3 - 


e 

c • 

o — 


0 — 0 * 


o o 

o o 

fc- w 

o o 


o ox 

o o«« 


X o o 


o — — — c 
— woo# 


o c 

C | 


o « 
o 


V o — 


— ^ o 


O W- woe 

0 — *« f o O — O w 

— w £^9UO£« 

— — 0*30— CXOO. 

1 O O — — ** — vj 

•“ f ? • xcoo 


o 

c 


CM 

*n 


a 

1 


o o 


— o 

• o 


x ac a 
— o 


3X00 
3 | a.— 
o o O X 
O X w o 
> O CL > 


— a— o 
, x x c a 

— O — o >s 

— © g +* 

• CO 

o o — a. a 

3 UJ O CL CL 


we 6- CO o o 

#wow%.o — — ec— o* 

o « o • c“«« S.C o o- O.B 2 S ° ° 2 - • 

-C 'J 'ocxxo coc-S .?- * !«oi!*.ocfoo QO !o'‘S**^ - X © # 

^il=l=illf IJsIIs jii : ?* : : H ff IslI t j l! « :i= 

lliliillj!lj!iflllliiill!«||g”|!|||||||||| 


*o 

© — 

• w - o 

• o o — 

** — — X 

• 0X0 


— — c 

<30X0 
c w o O 

w C C C 

o o 

X — 

3 c a a a. 

**# 3 3 31 
XXXI 
0— 0 0 0” 


— — o aaa 


wo — — ** 
# o o o o 

0.0 0 0 0 

xxx 

&«. 

3 O 

X o o o 
o K C C C 

— a o o o 
o — ** ** 

— X O O O - 

o—ooo 

O c w k. I. 
X Ul Iw Iw Iw I 



sis ••• •••« 

ZM0C#O#*CMO^O0#4 


• CM A 


> — lO®— *C4 — — 


• IO M 
• W^W 


x « o- 

OlOh^iOfO# 
— k »a ® in ® 


*o 

— rw < 

— — cs 


— o 

« CM *n fw 

*0 — ® — ® ® c* 
* CM CM *1 ® ® ® ® tf) 


4-120 







3«3$5*$S 


«o 


«-* r* i 

o«o« 


- NA NfNfNOC* A 

s :!sns:s . •sisisaaf.sf-s 




i«»N 

I © ^ © 


4-121 



o 


m - 


c 

A c 

k- O* — - 

• 3 

— o o • 

— k- — 40 - 

Q. •** C 

— • o • 

40 o »- c 

— O 9 
3»-£ • 

ft k. 40 • 


C • ft 


CL k- CL 
TO O O 


c 


• • 

• 


c- 



• 

1 C 

40 
40 • 
C 3 

• 

A 

c 

O 

C • k. 

— • o 

; -^ 
C X* k 

— c 
N O 
N — 

o • • 

c c « 

*0 

• • c O 

C 40 O O • 40 

— »-'O 0 • 

— X 40 — ■ _ 



• — » 

* U TO 
** O • 

40 

A 

•A C • 

— c o 

— O A A 

4 - c 

• V 

4# 4* 3 

c x *o 

C k» 
3 O 
O 4* 

• o o >s k. 

K 40 O • # 

o O C O 4# 

• 

c 

o 

— c 

• A 

3 >s C 

» u o 

> C Ok. 

C O A ** 3 

O 40 C 4* 

-* o + 

C 9 4- 
Ok. • • 

• • 
d • 

w • o — 

« o • 
• • 

40 — * c 

A — • CO 

«« t k o D • 
• o 5 o — w • 

• 

o» • 

o • 

40 Ok 

40 o 

o — 

C o 40 

A L. M 


• » • 

c — 9 

C • 0 • 

O C u 
£ C 

O 0£ v 
X 4- O 
9UO. 


_ o 

• “ V** • • 

oac T)h “D TJ 

• - o — 

• ~ — » o 

— «* e o 

C 3 « 

1 O C C T> ++ 40 

• o c e 

^ • o o c c c 

a X • • • 

«- o u o o» o» d(L 

o • • • </) 


o o» < 

O I 


• • p 


- o 


«• | 9« « O C Ov 

a S e c e t » • 

1 >s 

W O. 

fc. k. k w w 


m ~ 
o • 


9 • • 

■ o c 


• AC 


I- — § »- — 

• TO O • 3 A 


u c «- 

w O O 

• k- — C 

«« w «. • o 


o 

• 

k- o 

o • 


o • o 

• 3 


• 3 A 
«* w w 


— • A 
0 OX I 

O • o 9 
OK2Z 


Z'C O t I — ** • 3 o • N > > 

V • c • • k. 4* k. • m — C N 40 ■ l» 

' ** ^ 0^00000 x (% a 5 C C c a eft ^ d Z” 

^o ; QCk.oft i0 k.ae>5eio:8;a#“!.?L*2 

o 4 - • *. .. .. .. .. .. M 

* * * »••••••.!•••.; •••••• •••••••* 2.g Z w 

° aaaaik.X&& & a a* a " aT aT1T~TTT | A — 4- 40 

o • 5 aw • • 

I ftl 3 3 X f 

_ _ o 3+*+. 0 5 

«>.e 0 2222 2J22 222S sa » 3S *3*sss3"-»-a51! a 

•” w 9 9 ■ ■■■■■■■■■•■•■•■■• • ■■■ a o o « o ~ 

? ftO • 03A03A03A 

i. k, t e*« w w e*# w 

“ 3 — 3 *— 3 

-4«A-«#A-^ 

• k. # k. • 

3X33X33X 
k- O ** v O « *• 0 


— •«.! a. a a a a a a a a a a a a a. a a a " " " “ ° — 
•*• x« — ___.__rr a a a a a o • 

o o 


_ r“““« aaaaaaftaaaao • o 

-•©r~rrrrEEEEEEEEEEE~~E~~-~?-®*" a 5- 

»t °«iiii*« j sss *’’* 3 *’’***’’ — 

! | f f i ill l 1 1 1 ! f iiimiiliiiiiiii iiiiiiil i I i| ! 


»- x- «. *p 

o o o 


i-c~ 
o • o m _ 

— m — x» — 

t> o • o c 

— ■ *o 
#to — • o 
c me 

- to- • 

A C O A C 


• 033 


3 O C 
— C 
3 0 0 
c c 
c — # 

O A C 


A A A A * 

ms 

C C C«< 



• • M 

— • *0 tO *> 


4-122 





B 

X 


S 

• t- # 

• ft t > 
O ft 

+* 9 *o 

ax • 9 — 

S 9 ~ I o 

a • x« — 
o ft • xe 
X • o 

»* o — c 
9 0, W C 


IV 

k» 

o 


** ft 


• o 
ft «* 


o ^ ~ § 

M 90 xo i 


x w ft — 

9 ft *> 

— 9^ • 
• C ftft 

MCI 

X W k 
X O 9 • 

o ft X • 

X ft ft 3 


»D ftl 

w o w ** 

SSI- 

— ** c 
• ft* t- 


X • o a 
c o * « 
o o o u ac 
«« I- • o 

u . C k. • 

TMi • ao«« • 

~ c v ■ k • 

o- c** w • e 

** w ft Oft x 

•. u k w«» • ■ g 

** o e • — 


-o-o — w 


>s 

ft > 


** C O 
— X 3 
ft DftV 

c— e • 

ft X o x 

o ** o c 


o o 
«« 

o o 


ft ft 
ft ft 


ft X-** ft ft 


ss- * :s 


» • • ' 

• > > 
ft V — — _ 
««ft & w k > 


5 


! ; h 

99AS 

C ft V ft 
ft 3 C 

W 9 — «*•? 


k. ft ** 

Sii 


3 ft — 
&,** < C« 


ft ON> 

5.135 


ft X 

ft K W 

2882 


58 
m 9~zz 

ft ft ft ft 

e %. c 9 *o 

— ft — 


• - ft* ft X #x 

- f - J o^ - ji*" • • • • • ft* 6 
k«> - ft — or o«*£ c e c e • ft 
<o •«« — — — ftftft 

ft ft ft o x — 'ft — (re * - 

•.>ft2uft<#4*»-yy ft*; o 5 • » • 

k. WCKftV — ft ft — — — ftft ax • 
vvtioaftotc-e*'’* — > c x x c *» 

«. ft^««*. 9 I 9- ft ft O ft ft i ft +* 
ft w ft o 0— ft.ftuooo **ft »* 

«* a • a«* ft • ix ■ • ft 

>s ft X— ft ft ^ • +? 


ft ft 
ft ft 


ft ft o o X 

«#««««** ft 
0 9 0 0 «# 


. C C C 

ft ft • 

k fc» w i 

ft O ft 

ft ft ft 

sss 


ft c 

S C ft 

ft > 
* 9 — C 

V ft 4« 


C ft c 
ft ft ft 


x o ft o Oft 

ss3> .-2 


c 

o - 
— o 
+* c 
ft ft 
ft — 
ft ft 


ft c 

— o 
** — 
p ^ 
ft o 
ft ft 

- ft 


ft ft ~ 


ft ft 
ft ft 



o «« t k-ccccece*--. 

ktAlmzzzzzlz 




55 

ft ft 

w w 

aa< 


mill. 

• o ft 


c o 
• • — 

— ft ** 3 

— a e 9 

— xo e 
#«« <» o 


— ft 


o**- 


9*0 

c o 

ft w 


Sr 

ft • 


9 ft 
C ~ 

• a 
— ft 


»- © O w 


esc 
« o o 


c e x i 
e — u i 


• 

o 9 
ft ft k. 

c C ft o 
9 B 

9 9 O w 

c c ** o 

ft ft ft ^ 


k. e 
o 


** ft 
«•- 3 

a 


a<»-a 



SSS 3 o 8 = ’- 55 S 

§ 99 X 
U- < t 

Iklk 2 



•••••••••••Oft 



• m • 

ft d ft M 


• Nfr I 

N^ftftl 



109 

~ I I 



9 9 9 ft) 9 


9 9 

. * 9 9 9 

019 . . . 

ft • 9 9 9 9 9 
9 9 • • 9 • • — 


4-123 



o 

O L_ 

o o 


-* o x ~ 

5 > > • c 

o x o — o»— 

X U 3 a 

^ O O 

C — 3 Ok— 

O • * • o *© o 
CCO^*CCu 

— o # — a o — a 

He S 

• •oooooo*- 

OCk.ww^wwn 

o — ^ te ^ 

“ • I OOOOOOOOOXO 

• t S «««■■■- o«- O 

_ 2 4 s- '' £ c c c c c x%» o x 
® ^ OQOOOOO** f «,» 

g- ^cececcc o2 * 

'.Ji::::;;:-*,; 

rt 3 2*******0^l»3 
gJ^l'O’D'O-O-O-O-O^ C 3-2 

♦'■3jfbkm.ww W-2JS 

x jig:ss::::s K *i 

oai- 33 3;;;;’5i; 


x o 

*•* • 

— X fc* 


~ I «~ 


O I- • 

*- o — 
u ** 

OOP 

*- o c 
<•» — 

>1 

~ X C 

• *- 3 

^ • O 

a 

• p 

• • 

• e 

• -x . 


3 *■ o c — Ok. 

I? * • • • • 3 

U >**“ X 9 C« 

°o o pood 
c x c x • ’O • i. 

• o • o o • . • 

^ C — f | « ft 
» 000 «> 0 « 0 C 

•— • • 5 

*■ OT) > o «4^ 

^ — *" *“ • W | 

• CO « 

4 w *- o • - 
C t c cv o 

o O ^ O 3 1 . 

• O m %rnZ Z ° 

9«9t 3003 

a— a— o-— o 

K3M303VW 

• a* ao ol£ 

* * o o • 

» • • • Z _ 

C3M3X3KO*. 


o *< 


• - 
• o 


a 

• a 


o i» 
X o 
C- I o 
^ cx 

2 ° 1 
o c 
c — o 

t- • 

3 3 X 
X v o 


0 3 3 
l O O 


— — o 

— o » 

2 * 

» X 


o X C O K 
“ c o ** 

~ O m* — 

"O • 

*- V t c 
• V to • 

. ** • *o c *o 

>-o e • Ok 

| £ • a. c ° 

— CL • “ 


C 

O 


x • 


O C 


O 3 

3 © 

• — » k. > 

■O • • • 

»- c • 


c 


c 

?3t O • — 


O «« i. • « 
• ^ O k- C 


- ! 

m c 


X «# C ** 
O — 

* c 


— c e 

3 0 0 


C|K 

_ “ o o • • • 

'*■' <* • • t COM 


O 

3 X X M - 

•• e o o. 


2 Cl • w m © _ 
> o a3 i N 
*: o — c n 

O C C u o o © 
*- O O o > C C 

o o<*- • c c 

c k. o m - • • 
- - o o — ox o o 

“ ® ^ t c • 

X^«.OCJK3# 9 

a o-S 0 ^*"^ ** 

icccxtoic C^x C • k 

••••g*-OOOOC*-OUO 

Zl*Z~~ m Ji2e~ gg 

0 iiV 0u ^ # * f “>oo 

XXX k- — o k. XX fc- ** c • • 

t t • • » t • S 3 0 0 

• oSoSo 0 **!"***! 00 
|ao.a.*^"SS.S.Jco-f 

eitoQ.at|o« a- coo 


o o c • 


o • 
• 0 3 


C >, 

• — 

*» o 


o o — 

• • 

-C x o 

o o ■** 
— — o 


c m 
m 


O C II U g o 


a a a 

• • o 


9 o « 
c — o 


• X p 

o c o 
o o • 


c o* c 

— 3 X 

o o> 
** 3 

• o 

O • k. 

— o 

o • 

• «*. O 

O k o 
O 3 w 

• • k. 

X 3 


o o >s 


O O O 

a. 


o o x 
w k. o 


> > x x c 


c c » 
o o 


o o — 
o o o 

c c - 
o o o 
c c > 


C 3 3 © 


• • • 

• • m 


o — 
o 

— • 

• 3 X 

> *- o 


a. cl~ 


X K K 


ifrelpgl i8#I5|EgSigB8S“-- J * 

saiiilsaisssiftiSfttlii&siiBjiiijj 


■ tOjOoww: — 



K> 

> • <M ^ 

;• otom r 

2^222*° 40 « «- 
— « .fliS’lf ** NNlON# 


• M 


M 

2 2 *° 40 


► • ^-SS J< 


4-124 


N 





•- CM CSi K> 


A 9 

*- X 9 ® 

• • • • X ® 



«- tf) X X «- 


x 

X X 
• ^ 

— CM 



O • >t 


TJ 

O • 


o — - o +* 

b. — O • V 

a— 8 — 

P — o o e • — 

• o *# >- t •» ft 

*, *0 c O — 

c o ** e a o o 
a o • « - c a. - 

k«< b- ** - x** J 

30 U >b© O O O S 

41 O ft M L 3 

• Z m*-— w 90*0 • 

.OX #,,X 2 toi 3 
a~ 0 . o - b a« 
o I u ft o ft • 

«-~3 ft-occo 
o -•OOI'-OW 

*-k^Oft03-ft 
«« bft e o — A#* 
x a.«> b e — • t 

9 a- 9 c o c s 9 

« — c — — -ft o 

•caw«oepl- 

1— 03#^— w ®“ 
« «ft3 98 90 3 


«* 3 — - ft“ 

O • > • 

b fl — O • 

p «# O ft T) 

«# 3 — c 8 0 

p e b«< «« g ft 


*0 3 

o 8 
• 

£ O. 

-28 
o *© 

«* o 

o • • 

— * w O 

3 b. 

«> a u 

o • o 

• 

C b. o 

o a— 


8 

8^ 


^OOOOOOOOO 

*••••00000 


8 0 3 0 0 

• «. c b*.- 
U 0000 

— %. a. W o c *~ 

*. o so b ^ o 

— O a* 3 ft: <►* 

o x 8 *ft o 

a— ~ i i 0 

O > — *» o » 

- fl O C Xb 
«. *• c b o O b o 
0 — o — g clo — — ■- 

3 t I- 4 ' 4 ' • 0 

<Scwo< 55300 >b 

OOC 33 — ' O O 

S^ftft a #b 

• ScSSoO 30 

w •- «- - ■ o.j- • 5 » • 

1-22.0 

e 0.0 c c 

«# O b. — O 

OP — 

o o 

°82 


0 

0 — 

• -* 

• *0 0 


b — 

# — 

u 

0 
— * 

0 

ftfl b 

0 0 
a 0 


— 3 

- s 

a 

< 

a 

t> g 0 b 

a 0 - 2 

8 fOb « 0 

0 aao — — 0 

- f • t flb • 

X 
>s-* 
-* 9 

— >s 
+* — * 

3 X 

so 

for T 

0 ao— b. 
^ v g • fto ft 

— C 

O O 

0 k. 

x a 

m m 

0 
0 0 


>* 

♦» — ft 

>»— XO >s ~X 

** • +* C — » 9 ** 

- C~ O-X C 9 X 

300 * 00 — 'OC— ' 

c *9 C C 9 w • 9 

O O — O C — C 
■o — * 00*00 0-^0 
O _ w it b 


b. O O 
. - > b 

3 


O 

O O 
o b fl 

w 3 — 

o o c 

000 

«ri 

— a 9 


t- o o 


~ b- O 


a o 

_ _ o c o 

300«-0 —O — 

• b fl-- O X 

00000 — SOW — *k~ 

- 3 b. ^ — “ 


c 

T5 


3 fl C X O O b — • — o • — • * ** 

90 S .Ob^bflbt-oti*' • J* ; 

__.rz o o-ooo 89 - • 2 i\S 

oooa.>coo-*8^ o >'®5* T, 5. - .^ 
» *2 8 £■»«>>>*; 

b o 3 t^f L *faBl 4 beoc“ebbb-bCb-u-fl 3 i 

• a— c fto o- b f k 22 ’!S- »-- 4 -----fl^flfto 

• -ofl -oeir *ir >0t : . 1 T Z - b. - x ~ x 

3 K O — 


Sm»shmHs£umI 

3 * 




Ifooooooo 

g^flflOOOOft' 


b. < 

X 9 ft 
O — X- 


o o 

00 3 

3 3 a o 
O © S 
c c S ft 
o o as 

9 9 3 

o © — a 
«* — * o 

3 3 3 0 
O ««- O 


OO O- 

IftOfl 

ft | 3 o 

a 00 
o > > 

0 0—0 

ft o o e 
o — — 


o 

— 00 


-«< b b 
N O O O 


3 n C 3 

%- o — *•* 


*« -• too bJ co xft 0000 xft oaoftf > 

too>- 9 »oxooooxoooxoc # o- 

l i JS^.S5!-- # iS5S^ s ” s ” 5 2 ! 

.S(?22**SS>* .’5* 

IsiilSflEiii liflflflflflflflflflfll *ss 




JXU. 




ID K> I 
K> ID — 


i::!5fcJS:5!!!8iis8!::3j2llii::55:I5! 


• 9 r^lSSS 


99 tn ao 
9 ot 

9 9 »0 n fO CM 

if) M ® ® ® ® 


4-125 



o» 

C A 

• 

k- 9 
~ C 


o» • 
c c 

• • 

k- O 


X 

Ok 
— O 
. • — 
C %- 

"O 


c c 
O O 


o +* 

— c - 

o t 

o > c 

I • 

k- C > 
• O I 
*» — c 


o • 

A k- 
3 


C C 

» — o o 

o g 

JO 

o > a \a x x 


k- X • — 

#*««< o 
^ •- V- 

0 ft - k 

1 *0 c • • 

— —• — ** 

**•0x0 

3 X 4 • * 

c c • 

• • > 

O O 

00 o 

X X > 

-O £> O O 

i 

• • o 

c C • • Ok 
c c 

A A "Q 

k“ k- Ok Ok — 
3 3 C C O 
•'♦'•to 



• 

— 0 

• — 

A 

c 

0 — 

0 0 

O 

0 

0 — 

0 b. 

C 

b. 

k- 0 

a • 

— 

+* 

k*- 0 

• 0. 


c 

b. 

c 

b» 

0 — 

• 

— • 

• 

— 0 

b. 

0 — 


© o > > 


: c 
o 


• . 

3 

O *© X X~ « 

| C « o o 

O 3 3 

• o > > -© *o 

> C 

— C «*- +* O 

00000 
>—33 

*~ ^ O X 

k- o c c o o x<~ 
o — OOII^- 

~' 300 »-k__« 

0 • • • • C • 

“ c ££ C # 

3 — O O • x 

1 II T> 

3 Jtf b. b. M ^ 

o c • • U. U. « U. , 
o O A a O O -J Q - 


1 k- 0 3 


o . 
w 

"© • c • 

— b. I 

O. < 


0 • • 

D ■** Ok 

c o 

• 3 — 
> O — 
3 0 3 

• C -* 

o c 

1 — o 


O.C 


o c 
— o 


* o 
' c 

o • 

• ** Ok 

o 

I K ** 

o • 


3 • O 
• X— A 3 
k- •- — *» 

(3 3 0 O — 

O ■ X *< N-' ^ 

— •CO — 

*- • 3 • • O 

k- 0 £ O 

• — a. a c — 

• c X o o 

•'•OCk-O— — 
c 9 o O • A Ok 

X O «« b u c 

— W # «- «- o O 

O — Ok— c • 

£ > b O - b «. 

CSC 90 

• • & k.£X 0« 


X 3 

— o 


k *- c a— 


£ • o 
• O « o 


— c > 


b b •£ 

~^^* 00 >b. 

>31 «uiu o o 


o — 
o 


— AC 

* 8*5 


— X — 
— — • 

— o — c 

- • • • 

O C XI 
IX • 

C O C 

*0 0 Ok 

C «o 

O k- Okb 

• O A 

k- x b a 

• X£ O 

*0 o a>— 

X - ox 


k 3 e - 

O'*- — — O j 
— 0—1 
e c - b£ 

• • o *- 1 

• • k.x 

» ft • w 1 

♦'♦'XX • 1 

• • • C *0 ( 

A A X 0*0- 
C A O 4 

• • O — 

O O A — A - 
OO O < 
a & x 3 x : 

• • Ov o« 


*x>>5-H:<aj-^5dxlt2225^p? 
3>*«*a:2 0 588ll22;2£;^i£"<;5*§<g 



OQQOOKK 


)lvt 

► • 9 • 


1 1 ® — 

ujui 1 1 rs 

X X Ui bJ CM 

— A «0 CM — 

CMDfONtf 
A A 9 « 5 A 


_ . 

“•“•^••-™9AA99?CM®^. 


■ 9 

• O 

■ CM 4 


X 
I ( 
W 


■ ID 9 
>*> A I 


A A 


*NlA**0t(000A<A^00009^ 

4-126 


<2 




Output Listing 


4-127 




OUTPUT FOR MULTIPLE PUMPS AT FULL THRUST LEVEL 
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TOTAL ENGINE WEIGHT 13687. J4 



® ® < 

tO ® ® I 


► €> tf) 


O 

3 


o > 

CDhM 

x tr 

X O O 

UJ — i 
NU < 
CL 
Q 

— ui 
X -f • 
OOlL 
< UI 

o a 
X X 
O UJ K 
^ a ui 
■jxw 

U. Ui 3 


X ►- 
O X Q 
~ O O 
Ui • Qc 
* UJ 
* >- 
CO h- 
D O UJ 

O uj u 

Ui z < 
Z -i co 
< u. 
jujq: 

JQO 


O OC O 

v> ui < 

^ CO UJ 
X 3 a: 


*o 

04 

04 


x 

o 


•ONNOtCf 

• r-<M<NO*<N®® 

• K>fD®^-K>®® 

*0 ® ^ I** 

0* 04 04 


oc 

UJ 

M 

-J Q 
-M- UJ 
< Z => X 
OC 3 < U. O 

Ui O OC 

N-J«nDqq 
-— < CO CO OC OC 
JQDU(A<< 



< U. O CO OC Q. 

»- Ui UJ 

MOOOCNOl 
OC UJ Ui ~ Ui 

UiZZJOKUO 

ZDSDXhMN 

^ACDuotoai 


to PO 
«r> • 

— 04 
«o 

n 04 


:: 


X K ou 
OX JU 

— o ^ o 

Ui « O -I 
x uj 2 ~ 

5 * ->2 

Ot- UJ 

~ 2 3 X 

*- O U. o 

Z OC Ui UJ 

OD 12 

^03 I— « « 
u u 

CO CO 
CO CO 

ss 

o o 


o o 

-J J 

is 


o 

< 


o 


> 


o 

e 


o 

3 


o» 

to 

• 

• 

o» 

00 

K> 

04 ® ® — 

X ® ® ® ® ® a» 

*> 

• 

m 

\« ® ® ® co ® 




x«oo$t 

01 



O ® lO 04 40 GO 

— j m k) 


<r 

UJ 


CO 

z 


X Ui 

X 

< -j 


oc « rsi I o 

UQN t- ± 

►“ 0X0 

Ui 

UhZh2 

_l 

X ~ o UI 

<XU.2J 

I 

muouj 

ui 

O JQ 

> 

UJ QC < 

UJ -1 UJ UJ O 

_j 

ONfflO j 

< 


i- 

k o 3 h- < 

o 

CO Z Z to Q. 

►- 


CD 

3 

a 

Ui 

oc 

UJ 

> 


s 

I 

oc 

o 

u. 

oc 

Ui 

Q. 


o 

UJ 

CO 

o 

UJ 

oc 

UJ 

> 


X 

3 


< -» >- 
OC o Ui O 
Ui O Z 
3 <00 Ui 

H O Ui -3 ^ 

z < KUDO 

< >- < X 3 ~ 

-J * < -OU. 
u U «- LJH < u. 
UJ (O CO -I > Ui 
QLDQ.NZ • 

O oc *ni o: cl a. 
oc x a o 3 co co 
a. *- a z cd — -« 


4-136 



< 


8 

oc 

<L 


eot 

0 r* 


o 

Ui O 
to UI 

\<o 
-I ID 


I- UI 

<w 


1 % 

OuOC 

CL 

<ui 

SB 

*“ o 


to 

OC 

X 


3 

o 


s 

to 

s 


O 

0. 


28 

►" 0 
<0 0 

0 

8 

L. -I 
UI 

(0 > 

Ui UI 
_i —I 

O H 

UI (O UI 

X?j 
OKU 
(AX> 
I- o 

UI 

CEO* 

DUO 

$35 

8j uj S 
t *at 


CL K Z 

!Ls 

K< to 

n* 


a 

»— 

8 


* 


oc 

0 



UI 

N— ' 







OC « 

0 0 I 

0 1 

« 

00 

o 
o — 

0 0 

1 

0 1 

0 0 

Ui x 

01^ 1 

0 

0 f* 

o o 

0 — 

— 



Si 

IO 

«8 
a. u. 


x 

*- 

< 

tt 

ui 


N0 0 1 

0 0 * « 

* 0 Ml 

0 0 * t 


-1 

0 M 1 

000 — — M*0 

O 

0 0 


u. 

0 * I 

0**Kl^000 


0 

0 0 


< 

QC 


10 

(O 


a: 

a 


5 

— f 
UI 


< o 


to 

a. 


0 0 0 

0 

0 0 

0 0 M 

M 

M 0 

0 0 0 

0 

0 0 

0 


0 

* 




II 


0 

to 

UI 


0 0 0 
0 0 0 


UI 

0 

1 - 

0 O 

X2 J 

3 


• M00 — — — — 
^••lONNCSS 

0000Mr«0r^ 

— M M 0 0 0 


h h U. 

UI h- UI UI H- 

hJHU JhdCU 
ZhU Jh WwJ 

28*28 


0 

0 


— 0 

0 0 

0 

1# 

h* 0 

0 M 

0 * 

0 


3 ui ui 
O-i* 

UI (L UI UI o -I z o 

8||i;5ui|S“ 0 
£ 2 2££§*5S8 

h- HXUlhK 

*8252*8*22 UI 

<Q<<<MUIUIUJUJ u5s 
HOlslHKflCCfle Ohh 


28 
UI UI 


o 

UI 

o 


0 0 0 


0 0 M 0 — M 0 

-••0-00 
0 M M 0 
0 C4 


0 M 0 
0 — 0 


0M0 — *0M00 

000000 *00 
0 0 0 0 M * 

MM — 


£ 

O 



HZUIQC^ 


DtfjQZZ UI-IQ.0 

3ow«wwo»fcSg 

oOui< <^uiujp5 


28*32;:gii!8i 


4-137 


FLOWRATE SCHEDULE (LB/SEC) FOR STAGE #1 
GAS GENERATOR BLEED CYCLE 



• * • - 


CM 


CO 

SS2®' 0,0in ^« 

••• • • • • «o . 

• ® xo »- w _ 

*- <o 

HO 


9 

oo 


V) 

CL 


o § 

UJ tO (O 


v) a 

CL £*-s 

^oo 


UJUItO UJ ' — UJ 


tz a- o </> < j“r 


I 

K 

Ui 

IM 

O 


X 

O 


s S •*■ % ® x * o o — 

CL -JUI ft 2 K 

£3£58.- , "$£“!!L 

go o <K(t6»iT3d. 

w Jt§g^88wii!SS* 

psicssliim 

O.WWZZ<fl.>3j 


NND 1 

1 

1 M 


2*°® 


IO 


on ® 1 



® 

CM CM ® 

1 

1 

CM 

® 



HO KO 




*o m 



HJ 

♦ * 

m 


- 

m to 

® 


l-lh 
O O h- 
JOhU 
U.OUJ J 
to Z _l Z 

w CD Ui 3 
f dkm5 k 
o kuq: ui 

— < UJ y QC _j 2 

u-a>ujo<u.z 
£r3^r*<fflwa 

322-^8* *<, 


UJ < 

85 

wow 

J V) 
„ MOW 
: Nhoe 

> o £ 

> UJ z UJ 


CO o 


*5 


• o 
< 
DT 


W > 
O w 
< -J 


««• 

NN(0 

® ® rv 

m o 
« M 

m — 

CM ® 
® «- 
© «- 

® K 
® * 
® ♦ 

to »— UI 
to -J 

#30 (O 
© Dr >- < il 

67 

67 

64 

47 

17 

•n ho 

® HO 

« 

U- z o cl a 

^ H- 1-5 

>* or cl 





ori-Oi- 





< < h- U. -1 



UJ 


3 _ < < < 



QJ 


3<0CZO 



r\ 


S cl uj <o 3 


o o 

z 

u? 

to »- z u. 

UJ UJ — . 


o-ia 
CL O OH- 
KW CO * Z 

< -« Ul 
>- O <0 <J 

DC 


‘ * * ▼ • • • • • . 

•rt «cioS«io r 

N 


M 0 * 
M ® 

CM CM 


CM H) 


CM f 


(/) 

CL 

w 2^ 

0: ^to a 
p o < cl a.* , 

UJ CO CO •-» wo o 
UJ UJ CO to w wUJ 
CLOWCLOCWMs 



= IsE%!*sis|g= 
I rSg*Si c ‘5|ss 


V S h < £ m h n % 

g- SiwI^ 2 fS 25 {d. = 

alisiPsilSls k 

ilslslliliiil ? 


o» m « 

• • ® S ® ^ z ^ 

S2i«* s *-* 

QND *- 
CM CM 


CO 

a 


a: 
a 3 
w co ^ 
uj co < 

CL UJ •— i y— s 
V) Oc t/) (L 

a. cl z x-v 

O Z 

0 ^ ? woe 

1 QUOtrun^ 

KWOHWQUU 
wa w o to 2 l ~ h- — v 

P</)ql5wwowcd 
w to to or to <•• a -j 

<LKU.3> 

• z (O * - 

O O i 


® ® M 
® ® + 

® r-> cm 




S 2 

or i- 
U. < 
or 

• S2 ui 

*— < UJ OC 

■ to ~ o 
00 o to 

5 sr* 

CD 
or 
3 


_i u. or 
< w o. 


4-138 



» 9 IC» 

t M « O J< 


0 9 9 9 
M M Jt M 


m m c I . 

o o o z : 


9999999 
Jt J* It M J£ M J£ 


9 m 9 9 C 4 < 

9 K 9 9 
9^99 

«© -* r*. ci 

K) K> <N 

K> 

rO 


N 9 M'ttO 

9 * * * * ' 

♦ CM * lO •- 9 
900 -NiO 
« 9 — C* 

CM *- 


r^r^K>9999^«09 

0»K>I09999999 

h-tf>9K»9wr^«*-9 
CM K> CM 9 9 9 <0 
9 9 
9 
9 


* *■ r* 9 9 «- W 

♦ t©WB 

ION *“ ^ 

to 


9 9 9 9 9 K) 

* 9 9 9 I** OO 

tOMO ©»0 < 

r* CN CN *- 


9 . <M9 + - «0 9# 
9 «-)AN9iAn 
N 9 9 f* 9 
9 r* 

m CM M 


— 9 f 

— av — 


f 999 

9 9 9 <N 


11 e cl 


*9 9 

9 CM 



SJ1SS 

e 

9 * 9 CM 9 CM 

iDNf no 


M 9 9 « • • CM 
M 9 9 « » 
>N 

nmnwcM 

40 — * 9 

ioo<-^n 

CO 

9 


9 9 *■ 

K O 
UJ 


Ul (/> -j -I ui 

9 ^ AC 

N Ul 9 o a < 
(D Z 9 Ul t 
J m U U QC 9 
w socaow^s 
Ul Q£h 9 Uloa<J 
h ^3 9 Ul Ul -J P 
<ZK 9 CL Ul wZ 

S ~ hV)9LJ-Z 
wU. O Ul 91 < 
COOQCU 
-J Ul SkmUIWUI 
WhK 9 U. Z Z Z 
Ul Ul Ul mmmm 

SlSSSSiSS 

lollDflCtrthKK 


I 1= 

S SC 


Ul < 9 -I 
IVCOU. 
9 5 Ul 

1 ■ieigl 

-iuj»S%»-u!w 

SSJZgSwi 

3 -i Ul O X 

uj* 0 - 

yzixolo 

o ►- o o z z »- 


M UM>>U 
Ul ^9 Ul a 
H-< 9 Ul Ul p 

jpgall* 

lui<~UJUJX 

I-03UIKK0 

OOU.tfMQ.iO 


Ul« 9 
9UIZI- 
9 Ul X 


Ul ^ X >- 9 9 5 OJUj 

gtigsaS *sri“ 


9 ui x ~ 

Ul HO _ ^ 

jXXm 1 h 

N 5 Ul Ul X 3 9 

MM 9 Ui z> z> 

O 9 Ul htfXOK 

2S-JUI y uj h < I 


h- >- Ui 9 

X X 9 ►- 

o o z 

*Z*~UJ o 9 ui 

Ul Ul > -I ~ X 

99-1 O UJ < 

<tfu 9 9 

X Z > O ■— hOM 
“ “ Ul Z X UJ 


Isss^sssss 

U 9 < M XX 

X Ul Q X O UJ UJ ul 
uj -j Z -I < -j > > 

>N 0 << 0»-N-M 

y znohim tf mnjj 

jSoyotftxoyy 

NOZWH<hUJZOO 


1 J8^hX Ul Z X UJ 1 

lEig* li s 5i' 

IXZ J2M O ►“ < 

IOC w< xt-xui>. 
I hOOO < 9 O Z ' 


,z UJ Z3(l M Z j 

>-<ywh 

:<(L J<Q * x P Z < • 
;3hOOX OK9Uial< 


4-139 



I >»l 


to 


X 

o 


a: 

< 

UJ 


< 

h- 

o 


JtJtZZ-XJti 


cn r*. — m- 

N»OK)rf)®tn^ 

® onntoKi^© 

<N - *- rO CM 

*> *n <0 lO 


o * 
O • 9 


_ a. 1 • 

m m 

Z -* Z *0 

o o 

• 9 to 9 

to 9 

• 0t® 

— <o 

V ON ® 

m r** 

BMON 

9 K) 

(OMBCS 


to 


CM 




e e 

.o .o 



e 

•O 

c « e c e 

O O 42 O O O O 

ac 

o 

— o o» 

42 • • • 

— a • T 2 

c c 

CM 

r> ao m to to o» r> 

• • to 9 

• 9 

5427 

15687 

11296 

5 

6 

384 
13991 
11661 . 

99 V 9 

9 9 tO 9 

• O o 

• * 
to 

♦ 9 

— M- 
M- — 


a o >- 

-i — j -j 

X -J X 
(O UJ tO X 

r ^ 


3 to 

Sx 


V) 

o 


0 < X 9 

1 JO 


999l“99 

K£ 0 °*XI 

x x ~ *-» o p o 

OP^-KOXZ 
X~~<< 0:03 
Ld UJ UJ fie 0Z < < 
>-9 9 >-l3 

W h- K H 

MJWXXUUJUJ 

tozzppuzx 


Ld 

ac 

3 

to 

to 


< 

K 


x 2299 ; 

— uj ujS\qc i 

o ►“ o 


!2S2 Id ^ ^ 

: ui Ld 

► 

K<330<< 

-SSffSSS 


: — j _ » to to t- _j _j 


to ac 

z »- UJ 

iiu 

«IIWI 

§ >- o ~ o 

sff?? 

pill 


UI (O 


I 

I K 

UJ 
X t- 
K UJ 

23 

UJ ** 
I Q 

I UJ 

: z 

!22 


Q — J — J 
< < 

-> a : ac 

-J UJ Ld 

> 

o 


••!•••• 

i i 

••••£•• 

••!•••• 

i*y 

...... 

•••••• 

••••:*• 




UJ 







o 







< 







H 







to 







8 

u 







ac 







3 







O 

X 







u 

10 

P 







X 







z 







a: 







< 







* 







2 

I 

+ 

• 

to 

to b- 



_ j 

o 

9 

to 

to z 



_ j 

>- 

9 

UJ 

UJ •— 



o 



z to 

z < 



u 

to 


X to 

X K 






O Ld 

o ►- 



UJ 

•“ 

1 

~ z 

— ' to 



X 



X 

x z 



►“ 

1 

1 

»- o 

♦- s 




UJ 

< 

to* 

o 




to* 

X 

_J 




p 

o 

— < t- 

-J >- 

>- 



z 


< 

< DC 

UJ 



< 

UJ 

9 -J 

9 b- 

X 



ac 

p 

-J 

UJ 

o 




< 

-J < 

* 3 

< 



o 

t- 

< 9 

z o 

-J 



UJ 

to 

0 C 

< UJ 




o 


3 * 

1 - p 

z 



z 


t- z 

UI 



UJ 


o < 

o >- 

p 



s 


3 1 - 

CC ffi 

UI 





QZ 

< 

ac 



o 

o 

9 

b- b- 

9 O 





to U 

CC UJ 

z 

z 


Ld 

«0 

< 

o to 


p 


ac 

* 

to 

U. < 






z to 

Ld 

o 

to 




P z 

ac 

UJ 

UJ 



1 

~ o 
to — 

22 

ac 

UJ 

o 




UJ to 

p « 

»- 

UJ 




Q UJ 

to* 

z 

p 




o 

to UJ 

3 

< 



a 

t— 

2 u, 

“S 

UJ o 

°5 

to -1 

O 

¥ 

UJ 

1 - 

10 

-J 


CM 


* < 

V) 3 


< 


**■ 

X 

o o 

UJ 

UJ 

z 


* 

UJ 

3 

QC X 

to 

< 

to* 


a» 

Ld 

“1 

K Z 

to < 

? 

i 



_J 

-4 2 

b- 

a. 

J 


l 

fsi 

< 

ft 





fS| 

— z 

o »— 

to 

o 


ac 

i 

X — 

o u. 

< 

z 


o 

< a 

X < 

p 

Ui 


4-140 




Table 4-6. Sample Case No. 5 
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OUTPUT FOR MULTIPLE PUMPS AT FULL THRUST LEVEL 
LI0U, ° OXYGEN - LIQUID HYDROGEN 
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TOTAL INERT WEIGHT 25635.45 
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PRESSURE AND TEMPERATURE SCHEDULES FOR STAGE |1 
FOR ONE PUMP AT REDUCED THRUST LEVEL 60000. 
EXPANDER CYCLE 
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STORED PRESSURANT (AVE) 0.05 
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Table 4-7. Sample Case No. 6 
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Nuclear Thermal Vehicle 
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propellant summary for STAGE #1 
PROPELLANT IS LH2 
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ENGINE SIZE. WEIGHT,* PERFORMANCE SUMMARY FOR STAGE #1 
BLEED CYCLE 

CONVERGENT NOZZLE IS REGEN COOLED (MILLED SLOT CONSTRUCTION) 
NOZZLE IS REGEN COOLED (TUBE CONSTRUCTION) 

PROPELLANT IS LH2 
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the following is the regenerative coo 
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FLOWRATE SCHEDULE (LB/SEC) FOR STAGE #1 
BLEED CYCLE 
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ENGINE SIZE. WEIGHT.* PERFORMANCE SUMMARY FOR STAGE #1 
BLEED CYCLE 

(USER DEFINED GG) 

CONVERGENT NOZZLE IS REGEN COOLED (MILLED SLOT CONSTRUCTION) 
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OUTPUT FOR SINGLE PUMP AT REDUCED THRUST 



PRESSURE A NO TEMPERATURE SCHEDULES FOR STAGE #1 
FOR ONE PUMP AT REDUCED THRUST LEVEL 6M0e 
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5.0 MODEL VERIFICATION/COMPARISON 


NESS NTP engine system design, Sample Case No. 8 presented in Section 4.0, was 
compared to past preliminary engine system designs to support in verification of the models. Since 
no past detailed ENABLER-based NTP engine system designs are available that incorporate state- 
of-the-art engine system technologies, a comparison to similar, but not exact, engine system 
designs was undertaken. The 75,000 lbf thrust, 1000 psia chamber pressure, composite fueled, 
27(X) °K (4860 °R) chamber temperature, 500: 1 area ratio nozzle sample case was compared to a 
similar Rocketdyne NTP engine system design and a past ELES-NTP engine system design that 
are described in Refs. 2-3 and 5-1. The past ELES-NTP engine system design example did not 
incorporate an integrated ENABLER reactor system design, but included a reactor system design 
that only approximated in matching engine system cycle parameters. 

Tables 5-1 and 5-2 compare the NESS sample design case to similar Rocketdyne and past 
ELES-NTP engine system designs. Tables 5-1 and 5-2 compare key engine cycle parameters and 
major engine subsystem weights of the NESS sample case design to the Rocketdyne and ELES- 
NTP designs, respectively. One key observation is that the NESS design exhibits a delivered Isp 
of approximately 1% lower than that associated with the other designs. This is because the 
integrated NESS model more accurately calculates nozzle cooling losses. It was found that film 
cooling of the nozzle wall was required to keep its maximum wall temperature at or under the 
acceptable limit of 1460 °R. Table 5-3 shows the effect of wall temperature on engine system 
performance as predicted by NESS. The ELES-NTP did not properly model this effect. It is 
unknown if the Rocketdyne engine design properly represents this integrated design effect. The 
reduced Isp also increased the engine system flow rate slightly to offset this effect. 

The NESS program also more accurately models the pressure and temperature drops 
associated with cooling the nozzle and reactor system. This corresponds to the difference in the 
cycle pressures, temperatures, and turbopump operating parameters compared to the other 
referenced designs. 

It is believed that the NESS weights (and size) more accurately represents the reactor 
system because its model of the reactor is sized to take advantage of heat captured by the coolant 
before it enters the reactor. Likewise, the NESS integrated ENABLER reactor system module 
more accurately determines the reactor system weight and size for a given design point, when 
compared to past modeling methods, see Refs. 2-3 and 5-1. 
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Table 5- 1 Engine Cycle Parameter Comparison* 


Parameter 

Rocketdyne 

SAIC - ELES 
NTP 

SAIC NESS 

Total Flowrate (kg/s) 

36.7 

36.9 

37.27 

Pump Discharge Pres, (psia) 

1,544 

1,538.3 

2,298.3 

Turbine Flowrate, % Pump 

50 

50 

50 

Turbine Inlet Temp. (°K) 

555.6 

555.3 

622.3 

Turbine Inlet Pres, (psia) 

1,412 

1,416.8 

1,969.0 

Turbine Pressure Ratio 

1.25 

1.295 

1.739 

Reactor Inlet Pres, (psia) 

1,130 

1.255.4 

1,132.1 

Reactor Power, (MW) 

1,645 

- 

1,587 

Reactor Core Flowrate (kg/s) 

36.7 

36.9 

36.2 

Nozzle Chamber Temp (°K) 

2,700 

2,700 

2,700 

Nozzle Chamber Pres, (psia) 

1,000 

1,000 

1,000 

Nozzle Exit Diameter (m) 

4.15 

4.15 

4.22 

Nozzle Expansion Ratio 

500 

500 

500 

Specific Impulse-Vac (sec) 

923 

922.8 

912.9 

Pump Speed (rpm) 

37,500 

34,913 

40,583 


* Rocketdyne uses their Mark 25 type axial lurbopump (4 stages); SAIC ELES-NTP used a 
single-stage centrifugal pump; SAIC NESS. Sample Case No. 8, uses a 5-stage axial pump. 


Table 5-2. Engine Component Weight Comparison* 


Parameter 

Rocketdyne 

SAIC 

ELES-NTP 

SAIC NESS 

Specific Impulse - Vac (sec) 

923 

922.8 

912.9 

Reactor (kg) 

5,824 

5.823 

4.783 

Internal Shield (kg) 

— 

1,523 

1,108 

Nozzle Assembly (kg) 

440 

421 

535 

Turbopump Assembly (kg) 

304 

104 

221 

Nonnuclear Support Hardware (kg) 

- Lines, Values, Actuators, Instrumen- 
tation Thrust Structure 

1,815 

1,264 

1,493 


* Rocketdyne uses their Mark 25 type axial turbopump (4 stages); SAIC ELES-NTP used a 
single-stage centrifugal pump; SAIC NESS, Sample Case No. 8, uses a 5 -stage axial pump. 
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Table 5-3. Effect of Wall Temperature on Performance* 


Wall Temperature 

(°R) 

Barrier Temperature 

(°R) 

Isp 

(Sec.) 

Fuel Film Cooling 
Fraction 

1460 

1630 

912.9 

0.03 

1800 

2106 

915.9 

0.03 

2000 

2429 

917.5 

0.02 

2400 

2892 

919.4 

0.02 

2800 

3418 

921.2 

0.02 

3000 

3651 

921.9 

0.02 

3200 

3864 

922.4 

0.02 


* Core Temperature = 4860°R (2700° K) 


The ELES-NTP reactor system weight was approximated by reading off a reactor power 
versus weight graph that can have some inherent inconsistencies. The increase in the NESS 
weight for the TPA is due to the more stressing operating conditions under which the turbopumps 
must perform to meet the increased pumping requirements of the NESS design. Likewise, the 
NESS design comparison example includes axial turbopumps, which are not used in the ELES- 
NTP design example. The increase in the NESS design nozzle weight is attributed to a more 
accurate nozzle weight calculation that has been embedded in the NESS program. The ELES-N i P 
design approach only estimated nozzle weight by using multiple program runs to represent the 
various design portions of the nozzle. These results were then summed together which 
approximated the engine weight. NESS now calculates nozzle weight using exact geometric 
equations from which weight is determined. 

The nonnuclear support hardware weight is somewhat higher for the NESS design than the 
ELES-NTP design. The NESS design weight is believed to be more accurate than the ELES-NTP 
design weight because it uses true design calculations derived by TRW, see Ref. 1-1, during the 
past NERVA program effort that have been adjusted for today’s technologies, as discussed in 
Section 2.2.5. Additionally, the NESS nonnuclear support hardware weight calculations are more 
representative of an NTP engine system because they include options, such as those associated 
with a gimbal power supply, that can be a significant weight factor for long NTP engine bums and 
a weight allocation associated with a lower pressure cooldown propellant coolant feed leg. 3 he 
past ELES-NTP nonnuclear weight was estimated, based on a percentage of the reactor weight 
typical of the NERVA flight engine. This method has a greater degree of uncertainty. It is felt that 
the NESS program accurately models representative designs of solid core NERVA-type NTP 
engine systems to support preliminary design and mission studies. 
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6.0 CONCLUDING REMARKS 


The NESS preliminary design analysis program characterizes, in detail, complete near-term 
and next-generation solid core NERVA-type NTP engine system in terms of performance, weight 
and size, and key operating parameters for the overall system and its associated subsystem. The 
NESS program incorporates numerous state-of-the-art engine system technology design options 
and design features unique to NTP systems such as a multiple leg turbopump propellant feed 
system assembly and a low pressure cooldown propellant coolant feed system. The NESS 
program is easy to use and is flexible to address various NTP engine system design options 
efficiently. Though an initial validation effort, the NESS program is deemed accurate to support 
preliminary engine and vehicle system design and mission analysis efforts. 

Development of the NESS program is considered to be one of many key first steps required 
to support NTP development. Because of the modular nature of the NESS program, it has great 
potential for further upgrades in its design/technology option and analysis capabilities. 
Recommended future upgrade activities include: incorporation of other representative reactor 
system design modules such as for a particle bed and/or a pellet bed and/or cermet system; upgrade 
performance prediction correlations; include and upgrade materials option capability which 
considers radiation effects/compatibility; include a radiation heating model; integration of an 
efficient design optimization capability and perform more detailed analysis code verification. It is 
envisioned that NESS program could be a key design tool element when integrated into an 
advanced NTP engine system design workstation. 


6-1 




7.0 REFERENCES 


1-0 Re^ c i En ^rn 7 s ,rr/:r lation (NESS,; voiun " : 1 - pro * ram ««■» ou.de,- S ai C 

Report No. 265-070, NASA Contract No. NAS3-25809, December 1991. 

1 ' * £ le *? Uch ’ R d K ’ J R McDougall, R. B. Spencer, and K. R. Wener, "Volume IV- Detailed 
Technical Report, Nuclear Rocket Engine Analysis," TRW Report No. 8423-6008-RL000 
I KW Space Technology Laboratories, Redondo Beach, CA, March 1965. 

1-2 Taylor, C. E "Expanded Liquid Engine Simulation Computer Program - Technical 
si— ££££ R '- n N ° Oompany, 

13 uil° r ’i ? . E ' " Ex P anded Liquid Engine Simulation Computer Program - Programmer's 

“e ms'” 0 " NO ' ELES ' 1984 ' *"** ***— *—■ 

1-4 Livingstone. L M. and B. L. Pieree. "The ENABLER - Based on Proven NERVA 

, sol <on I T :eedlng ?^. lhe Elgh ' Symposlum m s P ace N “dear Power Systems, Vol. 
A pp. 598-620, January 1991. 

2 '‘ R ■ L P Dang ' ^ D E ColB - "Two-Dimensional Kinetie Reference 

Computer Program - Engtneenng and Programming Manual, "Software and Engineering 
Associates, Inc., Contract No. NAS8- 35931, Report No. SN63, April 1985 8 S 

2 ' 2 pXa,™ PCrf0rma " Ce Pred,C,, ° n a " d Cak “ la "°" **«-.■ CPIA 

2 3 oltr , Lun " and ^f s Mission Propulsion System Assessment Studies; Task 

Order NoJ. Nuclear Thermal Rocket (NTR) System Assessment Studies - Nuclear Engine 
ys em Weight Scaling Results." Inter-Office Memorandum from D. Pelaccio to M 
Stancati/H. Fetngold, Contract No. NAS3-25809, dated 9 April 1991. 

2 4 CPfiTth hr ’ T** L J quid P ro P ell “' 1 •"»«• Lland.ee Performance Evaluation Manual," 
CPIA Publication 178, September 1968. 

2 ' 5 Course SfS? uumA' ’ n ^^^'^teie m Technic,. M.,n„,l 

n TL m it a o £ 0 ' RJR> RockweU International Corporation, Rocketdyne Division 
Contract No. NAS 8 27980, 1 November 1984. y 0 ’ 

2-6 inie^tfonai r Main En( p^-" Sq a T ra n s i m aiinn teicni.rsgiu.kal Manual, Rockwed 
no ki TrmST??* 1 ' R «ketdyne Division. Contract No. NAS8-27980, DPD No it I 
DR No. LS-090-1, 1 September 1983. 

2-7 McCarty, R. D Hord, j. and Roder, H. M, Selected i kopertics of Hydrogen," U i 
Department of Commerce, National Bureau of Stands. KBS ! monograph 188, February 


7-1 



2-8 Liquid Rocket Engine Axial-Flow Turbonumpr, . NASA SP-8125, April 1978. 

2-9 Personal Communications with Dr. Stanley K. Borowski, NASA Lewis Research Center 
June 1992. 

2- 10 "Nuclear Rocket Program Terminal Report," Argonne National Laboratory, ANL-7236, June 

30, 1966, Published February, 1988. 

3- 1 Petrosky, L J "Scaling laws of Prismatic Fuel Nuclear Rocket Engines," Paper AIAA-91- 

2337. ALAA/SAE/ASME/ASEE 27th Joint Propulsion Conference, June 1991, available 
from the American Institute of Aeronautics and Astronautics, 370 L'Enfant Promenade 
S.W., Washington, D.C., 1991. 

. i * ' 

3-2 Petrosky, L. J. "Optimal Scaling of Pnsmatic Fuel Elements in Nuclear Rocket Engines," 
Paper AIAA-9 1-3506, AT AA/NASA/OAI Conference on Advanced SEI Technologies, 
September 1991 , availAbW from the American Institute of Aeronautics and Astronautics, 370 
L’Enfant Promenade S.W., Washington, D.C., 1991. 

3-3 Lyon, L. L. "Performance of (U,Zi)C- Graphic (Composited ind of (U,Zr)C (Carbide) Fuel 
Elements in the Nuclear Furnace 1 Test Reactor," I A-5398-MS, Los Alamos National 
Laboratory, 1973. 

3-4 Taub, J. M. "A Review of Fuel Element Development for Nuclear Rocket Engines," LA- 
5931, Los Alamos National Laboratory, 1975. 

3-5 McCarthy, J. R. and H. Wolf "Forced Convection Heat Transfer to Gaseous Hydrogen at 
High Heat Flux and High Pressure in a Smooth, Round, Electrically Heated Tank," ARS 
Journal p. 43, April 1991. 

3-6 Taylor, M. F. "A Method of Correlating Local and Average Friction Coefficients for Both 
Laminar and Turbulent Flow of Gases Through a Smooth Tube with Surface to Fluid Bulk 
Temperature Ratios from 0.35 to 7.35," International Heat and Mass Transfer Shortci 
Communications, August 1967. 

3-7 WANL "Radiation Heating Rates for the Typ$ C Components of the Preliminary Reference 
Ddaigs <939J723),” DRM-52852, Westinghouse Astronuclear Laboratory (currently 
Westirighotne AES), November 1970. 

5- 1 Pelabrio, ECC Scheil and J. Uvingstoi*v'7Jfciated Solid-Core Nuclear Thermal Propulsion 
Engine .Trades," AIAA Paper No. AIAA-91-3507. Presented at the AIAA/N AS A/OA I 
Conference on Advanced SEI Technologies, Cleveland, OH, 4-6 September 1991 . 


7-2 




REPORT DOCUMENTATION PAGE 


Form Approved 
OMBNo. 0704-0188 


Public .eporting ^Cen ,o, .N. carton „ SSS » 

^he^^nd marta.ung the data needed, and comping and reviewing J» >^^ 1 ^^ S( ^.^^a.e to. .Zmalion Ope-ations and Reports. mb Jedeison 

5g£asy^ papBmo,h ReduaK,n ^ <070<01881 wa5h,n9,on DC - 


“j AGENCY USE ONLY ( Leave blank) 2. REPORT DATE 

March TO 

4. TITLE AND SUBTITLE 

Nuclear Engine System Simulation (NESS): Version 2.0 

Program User’s Guide 

6. AUTHOR(S) 

Dennis G. Pclaccio, C'hrislinc M. Scheil, and Lyman Pclrosky 


3. REPORT TYPE AND DATES COVERED 

Final Contractor Report 

J I 5 . FUNDING NUMBERS 


WU-593-7 1 
N A S3 -25 809 


V. PERFORMING ORGANIZATION NAME(S) AND ADDRESS(ES) 


8 . PERFORMING ORGANIZATION 

REPORT NUMBER 


li-7623 


SPONSORING/MONITORING AGENCY NAMES(S) ANO ADDRESS(ES) 

National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, Ohio 44135-3191 


10. SPONSORING/MONITORING 
AGENCY REPORT NUMBER 


NASA CR- 19 1081 


11 S projecrMa T m^cr M °Im S cs T. Walton, Nuclear Propulsion Office. NASA Ixwis Research Center, (216) 433-4000. 


12a. DISTRIBUTION/AVAILABILITY STATEMENT 


1 12b. DISTRIBUTION CODE 


Unclassified - Unlimited 
Subject Category 20 


13. ABSTRACT (MtMlmum 200 words) 

•nm Program User’s Guide discusses the Nuclear Thermal Propulsion (NTP) engine system design features and eapah.l. lies modeled 
tins Program U« r • s Version 2.0 program (referred to as NP.SS throughout the remainder ol this 

7 ' ht ''f* 2 Klir^Ss ! JSL upgrade,, to include many new modeling capabilities no, available in the 

Keseai’cntcr December l,M. sec Ref. l-o. Ness’s new features inciudc: 

• An improved input formal 

• An advanced solid-core NIRVAlypc reactor system model (l.NAOl.l'.K II) 

• A hlccd-cyclc engine system option 

• An axial-turbopump design option 

• An automated pump-out turbopump assembly sizing option 

• An off-design gas generator engine cycle design option 

• Updated hydrogen properties 

• An improved output format 

• Personal computer operation capability 

Sample design cases are presented in the user s guide that demonstrate many of the new features ass.|ciatcd with this upgraded 
version of NliSS. as well as design modeling •caturcs associated with the original version «»l NLSS. d.scussed 


14. SUBJECT TERMS 

Nuclear propulsion; Computer program 

17 . SECURITY CLASSIFICATION 1,8. ^^ASS.RCAT'ON I,* 

OF REPORT OF THIS PAGE 

Unclassified Unclassified Unclassilici 

NSN 7540-01-280-5500 


13. NUMBER OF PAGES 

248 

16. PRICE CODE 

All 

| 20. LIMITATION OF ABSTRACT 


Standard Form 298 (Rev. 2 89) 
Prescribed by ANSI Std Z39-t8 



FOURTH CLASS MAIL 


National Aeronautics and 
Space Administration 

Lewis Research Center ADDRESS CORRECTION REQUESTED 

Cleveland, Ohio 44135 



Ml II I 



Posiage and Fees Paid 
National Aeronautics and 
Space Administration 
NASA 451 






